Static fluctuations of a thick ID interface in the 1+1 Directed Polymer formulation 
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Experimental realizations of a ID interface always exhibit a finite microscopic width ^ > 0; its 
existence is erased by thermal fluctuations at sufficiently high temperatures, but turns out to be a 
crucial ingredient for the description of the interface fluctuations below a characteristic temperature 
Tc(f). Exploiting the exact mapping between the static ID interface and a 1+1 Directed Polymer 
(DP) growing in a continuous space, we have studied analytically and numerically both the geomet- 
rical and the free-energy fluctuations of a DP at flnite temperature T, with a short-range elasticity 
and submitted to a quenched random-bond Gaussian disorder of finite correlation length ^. 

Centering our study on the two-point correlator of the derivative of the disorder free-energy R(t, y) 
as a function of growing 'time' t, we explore several analytical arguments at flnite f and propose 
a 'toymodel' in order to characterize the temperature-dependence of the DP endpoint fluctuations. 
Approximating the full correlator at small endpoint position y by R(t, y) ^ Dt ■ Ti-iy), we show on one 
hand that the amplitude Dod {T, ^) correctly describes the low- and high-temperature regimes along 
with the full crossover around T'c(C)- the other hand, we discuss the close connection between 
the temperature-dependent function TZ(y) at asymptotically large 'times' and the quenched disorder 
correlator. 

Finally we discuss the consequences of the low-temperature regime for two experimental realiza- 
tions of KPZ interfaces, namely the static and quasistatic behavior of magnetic domain walls and 
the high- velocity steady-state dynamics of interfaces in liquid crystals. 
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I. INTRODUCTION 

Effective one-dimensional (ID) interfaces can be spot- 
ted in various experimental contexts, encompassing do- 
main walls (DWs) in ferromagnetic or ferroic"*"^ thin 
films, fractures in brittle materials^ or paper^, contact 
line in wetting experiments^'^*^. The generic framework 
of the disordered elastic systems (DES)^^ has been proven 
to provide a quite successful modelling for such systems, 
describing them as point-like elastic strings living in a 
two-dimensional disordered energy landscape. The com- 
petition between the elasticity - the tendency to mini- 
mize their distortions - and the disorder - the inhomo- 
geneities of the underlying medium -, blurred by ther- 
mal fluctuations at finite temperature, accounts for the 
resulting metastability and the consequent glassy prop- 
erties observed in such systems. Moreover the value of 
the roughness exponent C, which characterizes the scaling 
properties of a self-affine manifold, is fully determined for 
a given DES once the dimensionality, the type of elastic- 
ity and of disorder are chosen, thus promoting the value 
of ^ to a reliable signature of the disorder universality 
class to whom a given system might belong. 

The specific case of a ID interface with a short- 
range elasticity and a random-bond (RB) quenched 
Gaussian disorder can actually be mapped on other 
statistical-physics models in the Kardar-Parisi-Zhang 
(KPZ) universality class""^^"^"*, including in particular the 
so-called '1+1 Directed Polymer' (DP) which has stimu- 
lated an increased activity lately, among both statistical 



physicists^^"^^ and mathematicians^^^^^. A large variety 
of results emphasizes the deep connection which exists 
between the descriptions of a wide range of systems up 
to random matrices^°'^^, such as the Burgers equation in 
hydrodynamics'^^, roughening phenomena and stochastic 
growth^^, last-passage percolation^^, dynamics of cold 
atoms'^'^, and vicious walkers^^'^^'^''. A shared feature 
between those related models is the well-known KPZ ex- 
ponent Ckpz = 2/3, which characterizes the exact scaling 
at asymptotically large lengthscales or 'times' assuming 
an uncorrelated disorder^'^'^^"'^'^. 

Although of interest regarding the whole KPZ class 
problems, there are two additional issues which turn out 
to be relevant for the study of experimental interfaces: on 
one hand, the characterization of the scaling properties 
at finite lengthscales, with possibly different regimes and 
crossover lengthscales regarding both the roughness ex- 
ponent and the amplitude of the geometrical fluctuations; 
on the other hand, the consequences of the interplay at 
finite temperature between thermal fluctuations and dis- 
order. However, in order to have then a complete realistic 
description, an additional physical ingredient must be in- 
cluded in the DES model: an experimental realization of 
interface always exhibits a finite microscopic width ^ > 0, 
which translates equivalently for a point-like interface 
into a finite disorder correlation length. Above a char- 
acteristic temperature Tc(^) > 0, thermal fluctuations 
simply erase the existence of such a microscopic width, 
whereas at sufficiently low temperature it becomes rele- 
vant even for the macroscopic properties of the interface. 
Those two temperature regimes can be hinted by simple 
scaling arguments^ ^, which are reflected in the two oppo- 
site Functional-Renormalization-Group (FRG) regimes 
of high-temperature'^^ versus zero-temperature flxed- 
point^^'"^^. Their connection has already been addressed 
analytically in a single computation in a Gaussian- 
Variational-Method (GVM) approximation^^''''*. Actu- 
ally its predictions for the low-temperature regime turned 
out to be potentially accessible and thus physically rele- 
vant for ferromagnetic DWs in ultra-thin films ^"■', which 
are believed to be the experimental realization of pre- 
cisely the ID DES considered here. 

In this paper, using the exact mapping between the 
static ID interface and a 1+1 DP growing in a contin- 
uous 2D space, we revisit several analytical arguments 
initially for an uncorrelated disorder (i.e. ^ = 0) in order 
to adapt and possibly generalize them to a correlated 
disorder (i.e. ^ > 0). In addition, we address the two 
issues of finite lengthscales and finite temperature nu- 
merically, by measuring both their geometrical and their 
free-energy fluctuations as a function of lengthscale or DP 
growing 'time', at finite temperature and in a spatially- 
correlated random potential. The first motivation is to 
obtain a numerical confirmation of the GVM predictions 
for the roughness of a static ID interface*^'"''*, at least 
qualitatively (if not quantitatively, because of arbitrary 
numerical corrective factors) : we actually successfully ob- 
serve the crossover to the low-temperature regime. The 
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second motivation, more general, is to test the validity 
and the coherence of the DP 'rounded' toymodel on the 
fluctuations of the underlying free-energy, as considered 
in Refs. 11 and 34, in order to correct some of the GVM 
artifacts, and to refine its assumptions so it would con- 
tain the essential features for the statics and hopefully be 
useful in order to understand the quasistatic dynamics of 
ID interfaces^^"^^ . 

An important preliminary remark concerning our nu- 
merical approach is that it is performed directly in the 
continuous limit of the 1+1 DP, and it uses an exact 
property of the model - the Statistical Tilt Symmetry 
(STS) - to focus on the effects of the disorder, disso- 
ciated from the pure thermal ones which hide them at 
small lengthscales. 

The plan of the paper is as follows. In Sec. II we define 
the full model of the static ID interface in the 1+1 DP 
formulation, along with the quantities of interest for the 
characterization of its geometrical and free-energy fluctu- 
ations at a given lengthscale r of the ID interface or 'time' 
t of the DP. Then in Sec. Ill we recall the exact prop- 
erties of the model at asymptotically large 'times' or in 
its 'linearized' version - obtained by neglecting the KPZ 
non-linearity ~, and use them to justify the construction 
of a DP 'toymodel'. Extensive scaling arguments are 
given in order to tackle the opposite low- versus high- 
temperature regimes and their connection, and the un- 
derlying scaling assumptions are actually made explicit 
using saddle-point arguments. In Sec. IV we perform 
an extensive numerical study of the free-energy and geo- 
metrical fluctuations of the full-model, whose results are 
interpreted from the point of view of our DP toymodel. 
After a synthetical discussion which yields an effective 
analytical prediction in Sec. V, we discuss in Sec. VI our 
results with respect to two experimental systems, namely 
the domain walls in ultrathin magnetic fllms and inter- 
faces in liquid crystals and we conclude in Sec. VII. 



II. 1+1 DIRECTED-POLYMER FORMULATION 
OF THE STATIC ID INTERFACE 

A. DES model of a ID interface 

We consider a ID interface, living in an infinite and 
continuous 2D space of respectively internal and trans- 
verse coordinates (z,a;) € M^. Restricting the model to 
the case where the interface has no bubbles nor over- 
hangs, each possible configuration is described by a uni- 
valued displacement field u{z) £ M with respect to a flat 
configuration defined by the z-axis (cf. Fig. 1 left). 

In the elastic limit of small distortions and for 
a short-range elasticity, the energetic cost of elas- 
tic distortions are given by the elastic Hamiltonian 
Hei [u] = ^ Jg^dz [Vzm(z)]^ wlth c the elastic constant. 

Assuming that we have a quenched disorder, account- 
ing typically for a weak collective pinning of the interface 
by many impurities, the microscopic disorder is described 



by a random potential V{z, x) with the corresponding en- 
ergy "Hdis [u, V] = dz V{z, u{z)). The disorder average 
O of an observable at fixed disorder Oy is then defined 
with respect to the probability distribution of the disor- 
der configurations V [V], which is assumed to be Gaus- 
sian, i.e. fully defined by its mean and its two-point 
disorder correlator: 



V(z,x) = 

^ ' (1) 

V{z, x)V{z', x')=D- 5{z - z') ■ R^{x ~ x') 

with D the strength of disorder, which quantifies the typ- 
ical amplitude of the random potential. The disorder 
should be statistically translational-invariant in space, 
and it is actually assumed to be uncorrelated along its 
internal direction z and correlated on a typical length 
^ > along its transverse direction x. Finally, we con- 
sider the specific case of a random-bond (RB) disorder, 
i.e. with a symmetric function R^{x) decreasing suffi- 
ciently fast to encode a short-range disorder and with 
the chosen normalization J^dx R^{x) = 1. 

At equilibrium and for a given disorder configura- 
tion, the statistical average over thermal fiuctuations 
{0)y is then defined with respect to the normalized 
Boltzmann weight Vy [u] oc e~^["'^l/-^ of Hamiltonian 
H [u, V] = Hc\ [u] + 'Hdis [u, V] (the Boltzmann constant 
is fixed once and for all at /cs = 1 so that the tem- 
perature has the dimensions of an energy). For a self- 
averaging disorder, a given observable must be averaged 
analytically first over thermal fluctuations and secondly 
over disorder (O) , recovering in particular a translational 
invariance in space. 

The choice of those different assumptions is explained 
in detail in Ref. 11. In order to compute the GVM rough- 
ness of such a static ID interface, R^{x) was chosen in 
Ref. 34 to be a normalized Gaussian function of variance 
2^^, encoding thus the typical width ^ as the single fea- 
ture of this correlator function (cf. Appendix A). 

B. Mapping of the ID interface on the 1+1 
Directed Polymer 

The characterization of the geometrical fluctuations of 
the static ID interface goes through the determination 
of the probability distribution function (PDF) of the rel- 
ative displacements V{A.u{r)) at a given lengthscale r, 
with Auz{r) = u{z + r) — u(z). The contribution of the 
combined PDF of thermal fluctuations Vy [u] and of dis- 
order V [V] can be disconnected by focusing directly on 
the fluctuations of segments of length r on the interface. 
As deflned in Fig. 1, such a segment can be mapped on 
the trajectory of a directed polymer starting from (0, 0) 
and growing in 'time' t in the 2D disordered energy land- 
scape described by the random potential V{t,y). The 
fluctuations of the DP end-point y{t) at a 'time' t, of PDF 
V(t,y), encode thus precisely the translational-invariant 
'P{Au{r)) at the lengthscale r. 
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FIG. 1. Left: ID interface configuration of displacement 
field u{z) with respect to the z-axis; definition of its relative 
displacement Auz{r) at a lengthscale r with respect to the 
internal coordinate z. Right: Focus on all the segments of 
the ID interface starting from (0,0) and ending at (t,y{t)); 
definition of the DP's end-point position y{t) after a grow- 
ing 'time' t. The translation table between those two rep- 
resentations is given by {z,x) -f-*- it' ,y') for the coordinates, 
u{z) ^ y{t') for the trajectory, V (Au(r)) ^ Vit^y) for the 
geometrical PDF and B{r) = {Au{rY) ^ {y{t)^} = B{t) for 
the roughness function. 



The energy of a segment of lengthscale r o ii, of tra- 
jectory y{t) connecting (0,0) to (ti,j/i), is given by the 
partial Hamiltonian: 



dt 



^-{dty{t)f + V{t,y{t)) (2) 



with the disorder distribution defined by (1). Integrating 
over the thermal fluctuations at fixed disorder V , the un- 
normalized Boltzmann weight of a DP ending at (ii,?/i) 
is then given by 

Wv{h,yi) = r^*'^"''l?y(t)e-«fe-^^*^l/^ (3) 

with the underlying four DES parameters {c, r',T, ^} . 
The connection between this continuous formulation of 
the DP, well-known among physicists, and its discretized 
version on a lattice with the solid-on-solid constrain has 
recently been properly established'^^. 

We restrict our study to the case where the polymer 
is attached in y = at initial 'time'. This choice corre- 
sponds to the so-called 'sharp-wedge' initial conditions'^ 
of the KPZ equation as opposed e.g. to the 'flat' ones 
where the initial position would be integrated upon. 



C. Geometrical and free-energy fluctuations 

We start with the definition of the relevant quanti- 
ties for the characterization of the geometrical and free- 



energy fluctuations, whose numerical study will be pre- 
sented later on in Sec. IV. 

With the following normalization at flxed 'time' t: 



Wv{t) 



dyWv{t,y) 



(4) 



we can define the PDF of the DP end-point, respectively 
at fixed disorder V and after the disorder average: 

Vv{t,y)^^^}^.V{t,y)^VV^) (5) 
Wv{t) 

and use them for the computation of averages for any 
observable O which depends on the sole DP end-point 
position y{t) (and not on its whole trajectory y{t'), with 
t' e [0,t]): 

{0[y{t)])v= dyO[y{t)]-Vv{t,y) (6) 



-oo 
oo 



{0[y{t)])^ dyO[y{t)]-Vit,y) (7) 

•J —OO 

and in particular the different moments of the PDF (5) : 

/oo 
dyyit)'' -Vvit^y) (8) 

/oo 
dyy{tf -Vit^y) (9) 
-oo 

Actually the PDF V{t,y) is known to be fairly Gaus- 
sian (though the study of its small non-Gaussian devia- 
tions encodes relevant physics'^'''^""'*''), in the sense that 



V{t,y) 



^-y^/(2B{t)) 



(10) 



with its main feature being summarized in its second mo- 
ment, namely the roughness function B{t) and its corre- 
sponding roughness exponent C,{t): 

B(0«,«.)4™ (U) 

a proper exponent ^ being defined only if a powerlaw 
can be identified on a certain range in t; this is typi- 
cally the case at large lengthscales, the beginning of this 
asymptotic regime defining the so-called 'Larkin length'"''* 
Lc. In absence of disorder, the DP is a Brownian ran- 
dom walk whose PDF Vt\,{t^ y) is then exactly a Gaussian 
function of thermal roughness Bth{t) = —. In presence 
of a short-range RB disorder, there is a crossover from 
this thermal roughness at small lengthscales to an asymp- 
totic roughness B^uit) ~ t'^^'^ in the 'random-manifold' 
(RM) regime of large lengthscales. A ID interface is thus 
a self-afhne manifold in these two lengthscales regimes, 
its geometrical fiuctuations being characterized by the 
scaling y{t)'^ - A{c, D,T,^)t'^< with the diffusive expo- 
nent Cth = ^ at sufficiently small lengthscales (extended 
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at all lengthscales in absence of disorder), and the su- 
perdiffusive exponent Crm = § asymptotically large 
lengthscales (obtained in Ref. 20, 28-30 assuming ^ = 0). 
Actually the existence of a finite width ^ > strongly 
modifies the scaling of the prefactor A{c, D, T, f ) and the 
roughness crossover, with in particular a whole interme- 
diate 'Larkin-modified' lengthscale regime for tempera- 
tures below TciC) = i^cOy^^ (cf. Ref. 11 and 34 or 
Appendix A). 

The PDF Vitju) and its roughness B{t) are precisely 
the quantities accessible experimentally via an analysis 
of a 'snapshot' of an interface configuration (defined as 
in Fig. 1); however only a single roughness regime has 
been observed up to now in ferromagnetic DWs, which 
are believed to be the prototype of our idealized ID inter- 
face (e.g. C = 0.69 ± 0.07 in Ref. 1). From an analytical 
point of view, additional information can be extracted 
from the fluctuations of the probability Vv{t, y) itself, or 
alternatively from its corresponding pseudo free-energy 
Fv{t,y) defined at fixed disorder by: 



Wv{t,v) ^ 



exp [-^Fv{t,y)] 
Fv{t,y)^Fv^a{t,y)+Fv{t,y) 
with the following conventions: 

Fv^o (t , y) = Fth (t , y ) + T log Wv^„ (t) 



^lh(<,y) 1 y 



cy 



Tt 



Wv^ait) = \/27rSth(i) - \ 2tt 



(12) 
(13) 

(14) 
(15) 

(16) 



The decomposition of (13) defines the disorder free- 
energy Fv{t,y), which fully encodes the integrated dis- 
order encountered by the DP up to a 'time' t. This con- 
tribution can be dissociated from the pure thermal free- 
energy Fth{t,y) because of the statistical tilt symmetry 
(STS) of the model, whose different incarnations are dis- 
cussed in Appendix B. Indeed, with the particular form 
of the short-range elasticity | (dty{t))^ in (2) and y be- 
ing a continuous variable, the effective disorder Fv{t^y) 
inherits the statistical translation-invariance of the mi- 
croscopic disorder 'P\V{t,y)] defined by (1). Its PDF 
at fixed 'time' (and similarly any functional of Fv{t,y)) 
thus satisfies: 



V[Fv{t,y + Y)] =V[Fv{t,y)] 



(17) 



In order to single out the y-dependent additive contribu- 
tion of Fv(t, y), we also define the random phase r]y{t, y) 
in a kind of 'random- field' formulation of the disorder 
free-energy: 



Vvit,y) = dyFv{t,y) 

1 / rv 
Fv{t,y) = - 



dy' Vv{t,y') +ctev{t) 



(18) 



where ctey(i) is a y-independent constant. Note that 
the STS implies that P [rjvit, y + Y)]=P [r?v(t, y)]. We 



may assume that the scaling of the distribution V [F, t] 
and V [ij^t] is in large part controlled by their two-point 
disorder correlators, on which we focus our interest: 



C(t, \yi - y2\) =[Fvit, yi) - Fv{t, y^)]" 

R{t,\yi-y2\)= vv{t,yi)vv{t,y2) 



(19) 
(20) 



which reflect explicitly the translation invariance, and are 
related by 



c{t,y) 



y py 

dyi / dy2-R{t,\yi-y2\) (21) 
Jo 



or alternatively by dyC{t, y) — 2R{t, y) using their parity. 

The connection between the fluctuating disorder free- 
energy Fvit,y) and the PDF V{t,y) with its moments 
{y{t)^) can formally be defined as: 



{ym = / VVV[V] 



Jdyt 



VFV [F,t] 



Jdyy 



^Fv{t,y)/T 



k p~iFti,it-y)+Fv{t.v))/T 



(22) 



Jdy( 



-{FtUt.y)+Fv{t,y))/T 



where only the y-dependent part of the pseudo free- 
energy, i.e. the information encoded in the random phase 
r]v(t,y) actually matters. However, even the roughness 
B{t) cannot be computed straightforwardly through the 
disorder average, this would require e.g. the introduction 
of replicae in a GVM framework'^^'^^'^^ (cf. Appendix A). 



D. Feynman-Kac evolution equations 

Since we work with a one-dimensional object (the ID 
interface, likewise the DP), explicit evolution equations 
can be written down for Wv{t, y), Fv{t, y), Fv{t, y) and 
r]v{t, yY^ ■ We can thus follow the evolution with contin- 
uous 'time' or lengthscale t of the effective PDF-related 
quantities at fixed disorder, and also of the mean values 

Fv{t,y) and 'nv{t,y)- 

However, no such closed equation for the correlators 
C{t, y) and R{t, y), the normahzed PDF V{t, y)), nor the 
roughness B{t) of course, are available. This limitation 
in the lengthscale renormalization of the disorder-average 
quantities is conceptually similar to the fact that the 
FRG flow equations'^^''^'^'^^''*^ of the disorder correlator 
R(^{x) (1) are truncated in a perturbative expansion in 
t = A — d (with the dimension d = 1 for the ID inter- 
face), an exact analytical description at all lengthscales 
remaining thus unsolved. 

At flxed microscopic disorder, in a continuous-time 
limit and at finite temperature, the so-called 'Feynman- 
Kac' formula^^"^^ for Wv{t, y), is a continuum stochastic 
heat equation with multiplicative noise^^'^^'^^''^'^: 



\Wv{t,y)] 









V{t,y) 



Wv{t,y) 



(23) 
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where the normalization Wv=o(t) is usually hidden in the 
functional integration jDyit) of (3). In order to clarify 
the normalization issues that arise due to the disorder, 
this last equation is rederived in Appendix C both in 
continuous and discretized 'time'. In absence of disorder, 
we recover the standard heat equation: 

dtVv^o{t,y) = ^^dlVv^oit,y) (24) 

whose solution at fixed 'time' is the thermal PDF 
'Pv=o(tTy)= 'Pth{t,y), i.e. a Gaussian function of zero 
mean and variance -Bth(i) = 

Moving in on the pseudo free-energy Fv(t,y) defined 
by (12) yields a KPZ equation with an additive noise^^'^''^: 



T 



dtFv{t,y) = -d'yFv{t,y) - - [dyFv{t,y)Y + V{t,y) 

(25) 

So the free-energy landscape seen by the DP end-point is 
a KPZ growing surface, whose disorder correlation length 
^ lies along the internal direction of the surface, whereas 
^ has been initially defined as a microscopic disorder cor- 
relation along the transverse direction of the ID interface 
or growing DP. 

As for the disorder free-energy Fy{t, y) (13), it evolves 
with a tilted KPZ equation: 



dtFv{t,y) ^^^dlFv{t,y) 
-ldyFy{t,y) 



1 

2c 



[dyFvit,y)y 



V{t,y) 



(26) 



with the new additive term stemming from 
-\[dyFt^{t,y)][dyFv{t,y)\ = -\i^v{t,y). ^ Apply- 
ing dy it yields finally the evolution equation of the 
random phase r]y{t,y) (18) itself: 



T 1 - 

dtilv{t,y) =Yj^lvv{t,y) - —dy[r]vit,y)Y 



dy 



Vvit,y) +dyV{t,y) 



(27) 



The disorder free-energy and its random phase encode 
all the information concerning the effects of disorder, so 
both Fv=o{t,y) and riv=o{t,y) are zero. In presence of 
disorder those quantities are moreover completely hidden 
at small 'times' by thermal fiuctuations: 



Fv{t,y) « Fv=_ 



,o{t,y)^ Fvit,y)^0, rjv{t,y)^0 

(28) 

whereas they completely dominate the large-lengthscales 
behavior {Fv{t,y) sa Fv{t,y) -I- cte(t)), the evolution 
equation (25) and (26) thus sharing the same statisti- 
cal steady-state at asymptotically large 'times'. Those 
disorder-induced quantities can be properly defined at all 
'times' (cf. Fig. 2), yielding in particular the following 
initial conditions: 



Vvit = 0,y)=Siy) 
Fv{t^O,y)=0 
Tjvit = 0,y)=0 



(29) 
(30) 
(31) 



Fthit,y) + Fvit,y) 




Fv{t,y) 



dyFv{t,y) 



dyFth{t,y) 




FIG. 2. (Color online) Free-energy landscape seen by 
the DP end-point as a function of 'time' or lengthscale t. 
Top: Graph of Ffn{t,y) ~\- Fv{t,y) (imposing Fv{t,y) = for 
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simplification): the thermal parabola Fth(t,2/) = ^ flattens 
and unveils the disorder fluctuations Fv(t,y), which sketches 
a KPZ surface in its steady state at asymptotically large 
'times'. Bottom: Alternative point of view with the graph of 
dyFv{t,y) = ^ + riv{t,y), where the random phase is pro- 
gressively revealed with increasing lengthscale. 



Considering at last the evolution of the mean values 



Fv{t,y) and rjv{t,y), at first the translation invariance 
by STS (17) trivially implies that Fv{t,y) = cte{t) (and 
{y{t)) — 0). Exchanging the disorder average and the 
partial derivatives dy^t, on the definition (18) and on (26) 
respectively, we obtain: 



r]v{t,y) = dyFv{t,y) = 



dtFv{t,y) 



2c 



2c 



(32) 

R[t,y = 0) (33) 



whereas (27) simply yields the consistency check 
dyR{t, y ~ 0) ~ Q. So the evolution of the mean disor- 
der free-energy is directly given by the sole two-point 
correlator of rjvit, y) in ?/ = at a given 'time' t. 
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As we will discuss at length in the next section, the 
behavior of R{t, y) at small \y\ corresponds to the curva- 
ture of the disorder free-energy correlator [t, y) around 
y = Q which fully determines the amplitude of the geo- 
metrical fluctuations characterized by the roughness pref- 
actor A(c, L),T, ^). R{t,y) has essentially a symmetri- 
cal peak centered at ?/ = 0, whose maximum is well- 
defined for a finite correlation length ^ but diverges in 
the limit ~^ (corresponding equivalently to the high- 
temperature regime). The connection between this reg- 
ularization at ^ > and the evolution of the peak main 
features, i.e. its typical width and amplitude Dt, will 
be the two ingredients of the DP toymodel constructed 
in the next section and tested numerically in Sec. IV. 

E. 'Time'-evolution equations for the two-point 
correlators R{t, y) and C{t, y) 

There are no closed equations for dtR{t,y) and 
dtC{t, y), but the combination of the Feynman-Kac equa- 
tions (26)-(27) with the Ito's formula yields nevertheless, 
as presented in details in Appendix D: 

dtRit, y) =^dlR{t, y)~\ {R{t, y) + dy [yi?(i, y)] } 
--^dyR3{t,y)~DR'^{y) 

(34) 

dtC{t, y)^- [C"(t, y) ~ C"{t, 0)] - f C'(i, y) 

-^^C3{t,y)-2D [R^{y) - R^m 

which would be closed but for the presence of the three- 
point correlators: 

R3it,y) = r]{t,y)^rj{t,0) (36) 

C3{t,y) = ~2[F{t,y)-Fit,0)] [F'(i,0)]' (37) 

Neglecting the non-linear KPZ term in the evolution 
equation (26) for Fv{t, y) is equivalent to neglecting those 
three-point contributions. The solution for the corre- 
sponding linearized correlator R^"^{t,y) is given in Ap- 
pendix E for a generic disorder correlator R^{y). It will 
be used in the two next sections regarding both analyti- 
cal arguments and numerical results, in order to discuss 
on one hand the expected qualitative behavior of the cor- 
relator R{t,y), and to precise on the other hand the role 
of the KPZ nonlinearity in the short- versus large-'times' 
and the low- versus high-T regimes. 

III. EXACT PROPERTIES AND 
CONSTRUCTION OF A DP TOY-MODEL 

In this section we first recall the analytical results for 
the asymptotically large 'times' DP fiuctuations, exact 



for an uncorrelated disorder (^ = 0), and examine their 
possible generalization for a correlated disorder (^ > 0). 
Armed with those considerations, we then define our DP 
toymodel by assuming that the asymptotic scaling in 
distribution of the disorder free-energy F{t, y) is mod- 
ified only at small displacements |y| ^ C where ^ > 0'^^ 
and remains a valid approximation even at small length- 
scales. Using scaling arguments, we explicit the relations 
between the DP toy-model effective parameters and the 
ID interface parameters {c, D,T,^} in the two limits of 
low- versus high-temperature, and extrapolate a contin- 
uous crossover between those two regimes via the tem- 
perature dependence of the GVM Larkin length'^**. We 
conclude this construction by sketching two saddle-point 
arguments for the roughness, which use either the large 
lengthscale t or the zero-temperature limit of as a con- 
trol parameter in order to accredit our asymptotic as- 
sumptions for a short-range correlated disorder (^ > 0) . 



A. Free-energy fluctuations at asymptotically large 
'times' 

At infinite 'time' and in an uncorrelated disorder 
(^ = 0), the distributions V [F] and V [77] are Gaussian 
and their two-points correlators are exactly known 

C^=o{oo, y) = Dao ■ \y\ , R^=o{oo, y) = Doo ■ S{y) (38) 

with Doo = cD/T. The Dirac (5-function of R encodes 
the infinite-'time' amnesia of the DP with respect to 
the remoteness of its initial condition t = 0, and the ab- 
solute value of C encodes the scale invariance of this 
steady state characterized by F{y) ~ y'^^'^ in distribution. 
This steady-state solution of the KPZ equations (25)-(26) 
for a (5-correlated V{t, y) actually yields the prediction 
Ckpz — 2/3 for the asymptotic roughness exponent^^ (as 
discussed later in Sec. HID). 

As shown in Appendix F 2, Gaussian distributions with 
those correlators are steady-state solutions of the Fokker- 
Planck (FP) equation for dtV [P, t] and dtV [77, i] at ^ = 
only at strictly infinite 'time ', imposing Doa — cD /T and 
the boundary conditions riv{t,y)\y^±ao — 0. The corre- 
lator of the random phase (20) then coincides with the 
transverse correlator of the microscopic disorder (1), up 
to the overall amplitude Doo- Note that the KPZ non- 
linear term — ^ [dyFv{t, y)] in (26) plays no role in the 
determination of this asymptotic amplitude, since its con- 
tribution disappears completely with the chosen bound- 
ary conditions. 

Actually at ^ = the distribution of our total free- 
energy Fy [t, y) itself, given by the KPZ equation with 
'narrow- wedge' initial condition, is exactly known at all 
'times' in terms of a Fredholm determinant with an 
Airy kernel^®'^^"^^; it is non-Gaussian and at asymptot- 
ically large 'times' it tends to the GUE Tracy-Widom 
distribution^^'^*, which eventually yields a Gaussian 
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function at infinite 'time'. Its second cumulant corre- 
sponds to our correlator C{t,y) (for Fv{t,y)) and is ex- 
actly known asymptotically as the correlator of an Airy 
process^^. It schematically adds saturation 'wings' to 
the absolute value (38) of C(^=o{t,y), which appear at 
~ (y(t)^), i.e. where the transverse displacement is 
defined by the roughness B{t) at this lengthscale t (see 
Ref. 46). 

In a correlated disorder (1) with ^ > 0, the distribu- 
tions V[F,t] and P[ri,t] are a priori not Gaussian but 
we can still focus on the two-point correlator R(t, y) prop- 
erties. We know in particular that its integral must be 
zero at all finite 'times'^^, with the possible exception of 
strictly infinite 'time'; this actually requires the existence 
of saturation 'wings' of the asymptotic C^=o{t,y), which 
are pushed to y ^ zLoo at t ^ oo. 

Inspired by the Fokker-Planck scheme at ^ = case, if 
we assume nevertheless a Gaussian Vq [F, t] of correlator 
R{t,y), we show in Appendix F3 that at infinite 'time' 
the correlator 




cD 



(39) 



defines a steady-state solution for the linearized FP equa- 
tion where the KPZ non-linear term \dyFy{t^y)\ 
has been neglected, imposing again _Doo = cD/T and 
the boundary condition f?v(ii 2/)la=±oo =0- Note that 
Ain t\ must consistently be Gaussian. Actually the 
full solution of this linearized correlator at finite 'time' is 
derived in Appendix E for a generic RB disorder correla- 
tor R^{y), yielding the following decomposition: 



cD 

~T 



R'"\t,y)^R^iy)-b^"\t,y) (40) 



FIG. 3. Top: Schematic graphs of the two-point correlators 
C{t,y) and R{t,y) (respectively in full and dashed curves) 
at fixed 'time' t > tsat, which suggests the generic decompo- 
sition (42)- (43); they display the two characteristic length- 
scales ^ and it in the y-direction; the dashed area below the 
central peak of R corresponds roughly to the saturation am- 
plitude Doo and translates into the slope of the intermedi- 
ate linear behavior of C (since dyC{t,y) = 2R{t,y) by (21)). 
Bottom: Corresponding PDF 'P{t,y) whose variance is the 
roughness B{t); the dashed area emphasizes the most proba- 
ble positions of the DP endpoint, which exclude the large y 
and thus the 'wings' of C or the negative bumps of R. 



with 



-.R 



(-1) 



y/^/B^Jt) 



dw vj^t 



-«)[»+j//V-Bth(«)]^(-l) / N 



dw w e L 'V ^'IR^ ' (vj) 

(41) 



where i?^ (y) denotes the primitive of the disorder cor- 
relator, all the rescaling is purely diffusive with as usual 
Bthit) — ^ and limj_j.oo 6''"(i, ?/) = so that the asymp- 
totic correlator (39) is indeed recovered. 

This result suggests the following generic decomposi- 
tion at finite 'time' for the exact correlator: 



R{t,y) = D^- n^{y) 



b+(t,y)+b-(t,y) 
2 



(42) 



J^dyTZ^iy) = 1 ^ J^dyb±{t,y) ^ l\ft (43) 



suited for the asymptotically large 'times' where we 
conjecture the function TZ^{y) to tend towards the 
microscopic-disorder transverse correlator R^ijj) at high 
temperature, with distinct corrections b±{t,y) corre- 
sponding to the 'wings' in C{t,y) and moving to large 
y with increasing lengthscale t as illustrated in Fig. 3 
with £t ^ \/B{t) as discussed in Ref. 46. 

Finally, the asymptotic amplitude Doo in (42) is pre- 
dicted to be cD/T in the limit ^ — ^ 0; however this 
prediction is non-physical in the limit T — >■ so it must 
break down for temperatures below Tc(^) = (^cl?)^/"^ as 
discussed in Ref. 1 1 . We will examine from now on the as- 
sumption that for T <Tc the decomposition (42) remains 
valid but with Dr^ = cD /Tc, justifying it with scaling 
and saddle-point arguments and testing it numerically. 
Since the full linearized solution is known explicitly, we 
know that such a saturation of the asymptotic amplitude 
can only arise from the KPZ non-linear contribution at 
'times' below the Larkin length, before the corrections 
h±{t,y) separate from the microscopic disorder correla- 
tor R^{y). 
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B. DP Toymodel 

We do not know exactly the distributions V [F, t\ 
and V [77, t], or even their correlators C(t^ y) and R{t, y), 
for a generic disorder transverse correlator (1). As we 
have just seen, neglecting the KPZ non-linearity in the 
Feynman-Kac equations, it is however possible to go be- 
yond (39) and actually compute at all 'times' the corre- 
lators (^'"^(t, y) and -R^™(t, y), starting from a generic RB 
correlator R^{y) (cf. (41)); they reconnect of course with 
the infinite lengthscale limit (39) but their correspond- 
ing corrections b±(t,y) encode a pure thermal scaling of 
the roughness, inherited from the small lengthscales and 
kept at all lengthscales^^. 

Taking an opposite point of view, we have considered a 
DP toymodel constructed from the asymptotically large 
'times' properties of the random-phase yyy. This con- 
struction is consistent at least down to the Larkin length 
Lc (defined after (11) as the lengthscale marking the be- 
ginning of the asymptotic powerlaw regime), since by 
definition of Lc the DP fiuctuations display a scale in- 
variance at 'times' t > L^, this construction can actually 
be pushed to smaller lengthscales, down to a saturation 
'time' fgat at high temperatures but as we will argue not 
at temperatures below Tc(^) (isat < Lc). 

First we assume that the effective disorder at fixed 
'time', Fy(t,y) and 'i]y{t,y), have Gaussian distributions 
accordingly to their linearized FP equation. So they are 
fully described by their two-points correlators (19) and 
(20) (translational-invariant by the STS (17)) and their 

mean values r]v{t,y) = and FY{t,y) = — ^_R(t, 0) by 
(32) and (33) (which play however no role in the compu- 
tation of statistical averages) . 

Secondly we assume that the random-phase correla- 
tor has a stable normalized function TZ, all the possible 
'time'-dependence being hidden in two effective parame- 
ters Dt and (,t- 



R{t,y) 



DfU^^iy), I dy-n^^iy) 



= 1 



(44) 



This form generalizes the decomposition (42) but neglect- 
ing the corrections b±{t,y). This is a self-consistent ap- 
proximation since those corrections and the correspond- 
ing 'wings' of C{t,y) appear at ~ B{t), and by defi- 
nition of the roughness it correspond to an improbable 
position of the DP end-point of decreasing weight 'P{t, y) 
as illustrated in Fig. 3. 

Finally we assume that the function TZ^{y) coincides 
with the transverse correlator R^{y) of the microscopic 
disorder, as for the linearized FP equation at infinite 
'time' (39). The effective width and amplitude Dt are 
kept generic though but under the asymptotic constraint 



D^{T,0 = fiT,0 



cD 



(45) 



where /(T, ^) is an interpolating parameter such that we 
recover the correct ^ = limit (38) with /(T, 0) = 1. 



In Ref. 34, we have obtained for this DP toymodel a 
set of GVM predictions for the roughness and the Larkin 
length, with R^{y) taken specifically as a Gaussian func- 
tion of variance 2^^ (cf. Appendix A). Those predictions 
are constructed centered on the fuU-RSB cutoff Uc{T,^) 
of equation (AlO). Assuming •Ct ~ C and Dt ~ Doo, the 
form (45) and the definition /(T, ^) = |mc(T, ^) yields a a 
self-consistent equation for the interpolating parameter: 



in 



T 



n 6 



(1-/), TciO^i^cDY'^ 



that connects monotonously the low- and 
temperature scaling of Doo at /(Tc, ^) ~ 0.94: 



(46) 
high- 



T^Tc-.f^ (47r)i/6 ^ ^ i?oo(0,0 - cD/Tc (47) 

T->Tc : / < 1 ^ D^iT, 0) ~ cD/T (48) 

and hence for the Larkin length (A9) and the asymptotic 
roughness (A7) beyond Lc'- 



asympt 



it) 



mo-' 



22/37^1/3 



2/3 



+4/3 



(49) 
(50) 



According to these GVM predictions, the amplitude of 
the geometrical fluctuations y{t)^ ^ A{c,D,T,Ot^'' at 
large lengthscales has a temperature dependence which 
is damped as sufficiently low T below Tc{£,) > 0, whereas 
the superdiffusing scaling C = f remains unchanged. 

For the DP toymodel, there is thus a physically deep 
connection between the fuU-RSB cutoff Uc, the ampli- 
tude of the random-phase correlator R{t,y) at small \y\, 
the Larkin length Lc and the amplitude of the roughness 
at large lengthscales, this last quantity being typically 
accessible in experiments. 



C. Scaling arguments 

Coming back to a full path-integral representation 
for the roughness and more generally for any average 
of observables depending exclusively on the DP end- 
point, we present thereafter scaling arguments such as 
sketched in Ref. 11 for the ID interface model defined 
in Sec. II A. We also refer to Ref. 60 for a previous ap- 
proach. Note that we systematically disregard the nu- 
merical prefactors in the whole section. Assuming that 
the random potential V scales in distribution consistently 
with its two-point transverse correlator R^{y) and that 
R^{ay) = a~^R^/a{y), the rescaling of the spatial coordi- 
nates and of the energy yields exactly for the roughness: 



B{r; c, D, T, - a' B{r/b; 1, 1, T/E, ^/a) 



(51) 



where B is the roughness function with adimensional pa- 
rameters, provided the scalings factors satisfy the two 
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relations involving the Flory exponent = 3/5: 

a = pl/3c-2/3;,)3/5 ^ ^ (^-1/5^2/5^)5/3 (53) 



E EE cayb = (cD^^^i/s ^ (acL»)i/3 



(53) 



Fixing one of the scaling factors to a characteristic 
scale of the model gives three possible choices, each suited 
for the description of a particular temperature regime 
(high-T, low-T and their connection), with the ad hoc as- 
sumptions on the scaling function B{f; 1, 1, T, ^). Firstly 
with respect to the temperature T: 



E ^T, b = 
B{r-c,D,T,0 



r*{T) E 

(T»T,) 



-eth(r)s^ (54) 

6h(r)25(^;l,l,l,0) (55) 



catches the high-T scalings if the function i3(f; 1, 1, 1, 0) 
is properly defined (cf. Sec. HID). Secondly with respect 
to the finite width or disorder correlation length ^: 



1/3 



h^r,{T,) (56) 
1,1, 0,l) (57) 



catches the low-T scalings if the function 
i5(f; 1,1, 0,1) is properly defined (cf. Sec. HID), 
with r*(Tc) = ^5/3c2/3£)-i/3^ Thirdly with respect to 
the Larkin length Lc(T, ^), defined as the beginning of 
the asymptotic 'random-manifold' regime"**: 

5-L,(T,e)=n(T)//(T,0' (58) 

« = ecff(r,0 = eth(T)//(T,0' (59) 

E^T/f{T,0 (60) 

B{r- c, D, T, = CsB{r/L,; 1, 1, /, e/^cff) (61) 

with /(T, ^) an interpolating function between the high- 
T and low-T regimes for both the Larkin length and its 
corresponding effective width: 

/(T,0) = 1, /(G,0-r/Te (62) 
Te(r,0) - n(T) , T,(0,e) = n(T,) = roiO (63) 

Ceff(r,o) = eth(T), ecff(o,0 = e (64) 

We can now focus on the roughness itself, and discuss 
the consequences of a powerlaw behavior at larges length- 
scales, which is known to be governed by the roughness 
exponent Crm'^' 



|. A behavior such as 



B: 



asympt 



(r;l,l,T,0 



(r»l) 



(65) 



without any other parameter dependence (this constrain 
can actually be taken as another definition of the Larkin 
length (58) for f = r/Lc) implies for the rescaling (58)- 
(60): 



-Basympt(?-;C, T',T,^) 



ic(T,0 



(66) 



An artifact of the GVM framework is that it predicts the 
Flory exponent of the model for the asymptotic rough- 
ness exponent; for the ID interface = | whereas for 



the DP toymodel Cp°^ = | = Crm'^' ■ So either the asymp- 
totic GVM exponent coincides with the Flory exponent 
of (52) and all the temperature dependence is cancelled 
in A(c, D, T, ^), or they do not and the scaling prediction 



A(c,z?,T,e) = 



=3/5T 



1/10 



4/3 



D 



f{T,0 



2/3 



(67) 



matches the GVM result for the DP toymodel (50) with 
(45). It is important to emphasize that the Flory ex- 
ponent I is imposed by the rescaling procedure of the 
full model of a ID interface, whereas the exact RM ex- 
ponent I is the true physical roughness exponent at 
large lengthscales and is predicted by assuming only 
that the scaling of the disorder free-energy is dominated 
by F{t,yf ^ D^\y\ as in (38) hence Cp"" = Crm'' (cf- 
Sec. HID). 

If we try boldly the rescaling 6 = r in order to catch 
the large lengthscales behavior, we obtain: 

b = r, a = (1)1/32-2/3^)3/5 ^ ^ ^ {cD^r)^/5 (68) 



B{r; c, D, T, (^) B(l; 1, 1, 0, 0) (69) 

that would predict the asymptotic roughness exponent 
(p^ = I if the function i?(f ; 1, 1, 0, 0) was properly de- 
fined, but this is not the case since the two limits T — > 
and ^ — > cannot be exchanged or taken simultaneously. 

The quantity /(T, ^) has been introduced here in order 
to interpolate between the two limits (54) and (56), in the 
only way compatible with the rescaling procedure (52)- 
(53). We argue however that f(T, ^) is the same parame- 
ter defined in (45) for Dqo in our DP toymodel. Actually 
all the scalings (54)- (57) are properly recovered in a GVM 
approximation of the Hamiltonian'^*, cf. (A1)-(A4), with 
the identification /(T, ^) = ^Vc{T,£^) that transforms the 
equation (A3) for the fuU-RSB cutoff Vc{T,£^) into 



167r 
~9~ 



T 



(1-/) 



(70) 



eA(c,D,T,^) 



So /(T, ^) turns out to be the key quantity for the connec- 
tion of our scaling arguments and the two sets of GVM 
predictions, centered either on the Hamiltonian or on the 
pseudo free-energy at a fixed lengthscale, both recalled 
in Appendix A. The numerical discrepancy between the 
equations (70) and (46) for /(T, ^) can be either reab- 
sorbed in the definition ~ |^ for the latest, or more 
safely attributed to the GVM approximation. 



D. Saddle-point arguments 

In conclusion of this section on analytical results, we 
present two saddle-point arguments which provide a con- 
trolled validation of our different assumptions at T > 
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and 5 > 0. Firstly we use ^ as a, large parameter at low 
temperature in order to argue the existence of a proper 
limit for B{f; 1, 1, 0, 1) in (57) (the high-temperature case 
(57) is already well controlled); this is not obvious in 
the usual conventions of mathematicians regarding the 
DP, see Appendix G. Secondly we revisit the original 
derivation of the exponent Crm'^' ~ f Huse, Henley 
and Fisher^^ from the point of view of our DP toymodel 
and using the lengthscale t as large parameter for the 
saddle point. 



1. Zero-temperature roughness of the ID interface 

The low-temperature limit in (57) can be made explicit 
coming back to the path-integral definition of the rough- 
ness and performing the rescaling (56) with t^{T) = 
as in (54): 



i?(ii;c,i?,r,o = e's( 1,1,^,1 



(71) 



2 Jyio)=o 



Vye T 

y(0)=0 



(72) 



where Vi{t,y{t)) = V{t,y{t))\o=i,^=i- In the path in- 
tegrals, the trajectories y{t) have a fixed starting point 
?/(0) = but a free endpoint y{ti). Since all temperature- 
dependence has been gathered in a single and large pref- 
actor the path integrals are dominated by a com- 
mon optimal trajectory y*{t), which, assuming that it 
exists, does not depend on temperature since it mini- 
mizes J^'^''''^''Ut[l{dty)^ + Vi{t,y{t))]. The saddle tra- 
jectory endpoint is then reached at some optimal end- 
point yl — Viii-jj^, Vi), common to the numerator and 
denominator and independent of T. Finally, one obtains 
from (72) that in (71) limr-^o i?(t;|5^; 1, 1, ^, 1) is fi- 
nite, being equal to j/i( f *r ) ' ^i)^- ^^'^ optimal 
path y* does exist and if its variance at fixed lengthscale 
ii is finite, the zero-temperature limit is well-defined. See 
Appendix G for a discussion on this last point. 



2. DP toymodel scaling argument, asymptotic roughness 
and Flory exponent 

The scaling arguments of the previous section, estab- 
lished on the full model of a ID interface, have of course 
their counterpart for our DP toymodel. The main as- 
sumption is that the large 'time' scaling of F{t, y) is 
governed by its infinite- 'time' correlator C{t, y) = Doo\y\ 
(38) with the amplitude being essentially a constant 



Dt ~ -Doo (denoted thereafter simply by D) and simi- 
larly ~ ^. This ensures that upon the change of vari- 
able y — ay and t — bt, the following free-energy is equal 
in distribution to 

2 -2 

Fth (i , 2/) + F(i , 2/) = y ^ + a 5 5 5 Fi (i , y) (73) 

where Fi{t,y) = F{t,y)\g^^ The argument of 

Ref. 29 can then be summarized as follows: the free- 
energy and roughness fluctuation exponents x E^nd C are 
respectively defined as F{y) ~ h^F{y) and y ~ 6^ (which 
amounts to take a ^ 6^). The fact that in distribu- 
tion F{t,y) ^ a^F{y) implies X = 5C while equating the 
thermal and disorder contributions in (73) leads yields 
X — 2C — 1. These two equations fully determine the val- 
ues of the exponents: X = | and C = §• Taking care of 
the prefactors of those powerlaws, we define the following 
rescaling procedure: 

Bit; c, D, T, i) = a" B{t/b; 1, 1, T/E, $/«) (74) 

a - (^/c2)l/362/3 ^ 5 ^ ^5-1/2^3/2 (75) 

E = ca^/b = {&b/c)^l^ = b-^l^a^l^ (76) 

where B is the roughness function with adimensional pa- 
rameters, if the scalings factors satisfy the two relations 
involving the Flory exponent (^^^ = 2/3. To understand 
how this power counting can describe correctly the large 
'time' asymptotics, we chose the rescaling equivalent to 
(68): 



b = t, a = {D/c^)l ti , E^ (DH/c) 



1/3 



(77) 



which implies from the definition of the roughness 
B{t;c,D,T,i)m (22): 



Bit;c,D,T,i) 



3 4 
t3 



/ dy y2 exp 1 - ^ 


. C 


1 
3 


f + A(1,27) 


1} 


/dyexp{-i[^i]^" 


f + A(l,y)" 


} 



(78) 



where the overline denotes the average over the random 
Fi . The advantage of our specific choice of the rescaling 
parameters a and b is that the 'time'-dependence of the 
exponentials in (78) is then gathered in a single prefactor 
^3. For each fixed Fi, one may thus evaluate the inte- 
grals in y through the saddle point method in the large 
t limit. The integrals at the numerator and denominator 
of (78) are dominated by the same y*[Fi] which mini- 
mizes ^ + Fi{l,y), ensuring that y*[Fi] is independent 
oft. We read from (78) that 



i?asympt(i; C,i?,r,0 = (^[^1])' ' (i^/c^)^ (79) 

i.e. the roughness exponent is Crm'^* = §• However all 
this construction breaks down at the very last when 
the scaling F{t,y)^ ~ B)\y\ ceases to be valid, at small 
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\y\ < when the scahng factor a{t) matches with 

the effective width f . This yields an alternative definition 
of the Larkin 'time' as a(io) = | or io = {c^£,^/Df^- 
Coming from the large lengthscales, this asymptotic scal- 
ing breaks earlier due to thermal fluctuations, at the 
Larkin 'time' tc > to. Identifying tc and Lc{(,,T), gen- 
eralizing ^ « ^ to ^es{T,^) of (63) and using finally 
D = f{T,£_)^ of (45), we recover consistently with (49) 
and (62) for the Larkin 'time': 

/ 2^3 \ 1/2 rp5 

a(t,)^^eff^tc= (^-g^j /(r,0-' (80) 

which has as a lower bound its low-temperature limit 

to^e^'cy^D-'/'^niT^iO) (81) 

The large- 'time' limit makes the scaling assumption 
Fi{l,y)^ ~ \y\ even more reliable, and the saddle point 
can be properly taken in this limit, yielding the Flory 
exponent of the DP toymodel Cf°^ = |- This was not 
the case for the ID interface in (69). Indeed, upon the 
rescalings (68) we obtain in a path-integral representa- 
tion: 

S(ti;c,AT,C) 

= S I; 1,1, r, r^2 s] 82 

_ "i'qi 5 Jy{0)=0 

A(o)=o 

(83) 

with a = (Di/3c-2/3^^)3/5 ^ ]^ ^ (ci:>2^i)i/5. The large 
tl asymptotics cannot be taken directly from this ex- 
pression since it is not in a saddle form and all scales are 
intertwined, contrarily to the study of the free-energy it- 
self which corresponds to scales integrated up to 'time' 



study of the geometrical and free-energy fluctuations of 
the 1+1 DP, and consequently of the static ID interface, 
directly in their continuum formulation. This approach 
uses an exact property of the model - the statistical tilt 
symmetry (17) (cf. Appendix B) - to focus on the effects 
of the disorder, dissociated from the pure thermal ones. 

Among the numerical procedures previously used to 
tackle this problem, we can mention firstly the DP un- 
der the solid-on-solid (SOS) constraint^^''*®''^^ where the 
polymer lives on a discretized lattice; secondly the semi- 
continuous ID interface, discretized along its internal di- 
mension but with each point living in a continuous ID 
splined random potential^^; and thirdly the continuous 
DP that we will present thereafter. Those approaches 
are of course complementary, especially for the investiga- 
tion of the large- versus small-lengthscales and high- ver- 
sus low-temperature properties, if a suitable translation 
from the specific numerical parameters to the physical 
ones {c, D,T, ^} is provided'^^'*^. 

Our procedure is actually based on the continuous ana- 
logue of the transfer-matrix method of a DP on a lattice 
with the SOS constraint, in the sense that it is performed 
after integration over thermal fluctuations since it fol- 
lows the partial partition function with the lengthscale 
(namely the Feynman-Kac equations (25)- (27)), for many 
individual disorder conflgurations. The continuous-limit 
formulation has two advantages: first the discretization 
issue in numerics is pushed back to a problem of numer- 
ical integration of partial differential equations, so even 
'small' lengthscales can be studied without discretization 
artifacts; secondly in the continuous limit the numerical 
parameters are directly the physical ones of the analytical 
model. 

Thereafter we will first present our numerical pro- 
cedure, and then discuss the different numerical re- 
sults obtained, which are respectively the 'time'- 
evolution of the correlator R{t, y) and its parame- 
ters {S,t,Dt\ at fixed temperature, the temperature- 
dependence of the asymptotic effective correlator Rsax ill) , 
the 'time'- and temperature-dependence of Fy(t^y), and 
the temperature-dependence of the roughness B{t). All 
the numerical parameters are gathered in Appendix H. 



IV. NUMERICAL STUDY OF THE 
CONTINUOUS 1+1 DP 

In the previous section we have considered several an- 
alytical arguments in the context of a finite disorder cor- 
relation length ^ > and combined them to propose a 
DP 'toymodel' which assumes essentially that the main 
physical features of the correlator R{t^ y) are gathered in 
ft, Df and the shape TZ{y). In this section devoted to 
a numerical study of the continuous 1+1 DP, we aim at 
testing those analytical results and the validity of our DP 
toymodel. 

The Feynman-Kac equations (26)-(27) for dtFv(t,y) 
and dtriv{t,y) provide the starting point for a numerical 



A. Numerical recipe 

Computing Fy{t,y) for individual disorder configura- 
tions V(t,y) up to a maximal 'time' tm, we can mea- 
sure directly the geometrical and free-energy fiuctuations 
stemming at thermodynamic equilibrium from thermal 
fiuctuations exploring a given random potential. Then 
averaging these fiuctuations over many disorder configu- 
rations we have access to the quantities of interest defined 
in Sec. II C and discussed in the previous sections from 
an analytical point of view. 

a. Finite box and microscopic discretized grid. The 
Feynman-Kac evolution equations (26)-(27) are defined 
in a continuous limit with t > and y CzR. In 
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numerics we work necessarily on a microscopic dis- 
cretized grid in both (t, y) variables and in a finite 
box (t,y) e [0,tm] X [— ym,?/m] without periodic bound- 
ary conditions. However, in experimental realizations 
of ID interfaces we also have some microscopic cut-off 
in lengthscales, ultimately the crystal parameter, and 
a macroscopic cut-off due to scarcer statistics at larger 
lengthscales. 

b. Generation of an individual disorder configura- 
tion. In order to generate individual disorder con- 
figurations for a given set of DES parameters, 
we first define a two-dimensional grid of spacing 

^^^f^,^f'i^={t^/Lt^y^/Ly} {d. Fig. 4). On each 

point of the grid we pick a set of random numbers ac- 
cording to a normal distribution of variance D^"*^ (the 
'strength of disorder') and the smooth random poten- 
tial V{t, y) is obtained by interpolating between the grid 
points with a 2D cubic spline as illustrated in Fig. 22 
{top). As detailed in Appendix I, its corresponding two- 
point correlator at fixed 'time' is given by the translation- 
invariant: 

V{t,y)V{t,0) = D"^ ■ i?cubics (84) 
D'=^ ^ D- i?g"^''=S(o) = £,s"'i^s"'^ (85) 

with the function ijCubicS ^^^^^ expHcit in (I2)-(I3). The 
values of those parameters are fixed for all our numerical 
computations to 

= 1 , ^g"d ^ 2 , L»g"d =4 ^ D = ^ 8 (86) 

Note that the physical parameter used in analytical ar- 
guments is D and not D^^, but they coincide with this 
particular choice ^f'"^. 

c. Numerical integration at fixed disorder. We have 
chosen to follow the evolution of the disorder free-energy 
Fv{t,y) (26) because among its counterparts r]v{t,y) 
is too noisy, Fv{t,y) includes Fth{t,y) that hides the 
disorder-induced fiuctuations of Fv{t,y) except at large 
'times', and the exponential in Wv{t,y) oc e^^^'^'^'"'^^^'^ 
suppresses the numerical resolution. Moreover, by sub- 
tracting the exact Fth(t,y) for y e M directly in the evo- 
lution equation, we get rid of the finite-box artifacts that 
would have arisen in the pure elastic contribution. The 
numerical integration of the differential equation (26) 
was performed using a numerical algorithm included in 
Mathematical'^, which adapts the numerical discretiza- 
tion in y at each 'time'-step in order to minimize the nu- 
merical errors. The main limitation was that the larger 
the maximum 'time' tm, the longer the computation time. 
Moreover, the lower the temperature the more the numer- 
ical solution of the Feynman-Kac equations is sensitive to 
the spatial variations of the random potential, thus dic- 
tating a smaller grid discretization in order to minimize 
the numerical error and increasing considerably the com- 
putation time. The number of disorder configurations 
that have been considered for a given set of DES param- 
eters is thus a compromise between the convergence of 



the disorder average and a reasonable computation time 
(as summarized in Appendix H). 

d. Initial condition. At short 'times' the free-energy 
is dominated by its pure thermal contribution (28) which 
diverges at t = 0, so we take as initial condition at 
to = 0.1 the thermal distribution Vito, y) = 'Pth(ioi y) as 
defined by (24), i.e. Fv{to,y) = 0. We assume that at 
large 'times' the DP has completely forgotten this ini- 
tial condition, but at short 'times' above its behavior 
should carry some artifact for the disorder-induced quan- 
tities Fy and rjv. 

e. Boundary conditions. We imposed at each 'time' 
that r]v{t,±ym) = i-e. (i, ±?/m) = 0. This is 
equivalent to the statement that for \y\ > ym, we have 
Wvi^, y) ~ Wv=oit, y) and thus the normalized DP end- 
point probability is 7^(t, y) ~ Vthit, y) ~ 0. This bound- 
ary conditio n can be physically correct only for 'times' 
such that yjB{t) < y^, else the DP 'senses' inevitably 
the boundaries of the finite box. This choice of bound- 
ary condition implies for the normalization Wv{t) (4) 
that the contribution for y outside [—ym,ym\ is exactly 
known analytically: 

\wr'\t) = j'j dy ■ e-(^-(*-^)+^Wt,.))/T + ,t,(^^ 
cte(t, y^) = y dy e-^-(*'^)/^ , F,^{t, y) = ^ 

(87) 

which yields Vvitjy) as defined by (5) (cf. Fig. 4). 

/. Recorded data. We have actually recorded, on 
a microscopic grid linear in y: on one hand the 
disorder free-energy Fy{t,y) and the random-phase 
rjv{t,y) = dyFv{t,y), with their mean values Fy(i, 0) 
and r]v{t,0) and their respective two-point correla- 
tors C{t,y) and R{t,y); and on the other hand the 
PDFs Vvitjy) and ^^{t^y), their corresponding moments 
{y{t)^) and in particular the roughness B{t) = {y{tY) 
with the roughness exponent C,{t)- We have recorded 
those quantities using a microscopic grid linear in t, and 
in parallel for the roughness-related quantities we have 
used a microscopic grid logarithmic in t in anticipation 
of powerlaws determination. Note that the even parity of 
the correlators has been explicitly used to increase their 
statistics, so the measured C{t, y) and R{t, y) are sym- 
metric by construction: 

R{t,y) = ^[T]{t,y)ii{t,Q)+j^{t,-y)7j{t,0)] 

C{t, y) = \ [F{t, y) - F{t, 0)] ' + [F{t, -y) - F{t, 0)] ' 

g. Fitting functions for R(t,y). The correlator 
R{t, y) is not known exactly at finite 'time', and 
we have thus only postulated its generic form 
i?(t, y) « A • 7^|^ (y) in (42)-(43). Following the DP 
toymodel assumption (44), we have systematically ex- 
tracted the typical width and amplitude Dt for three 
different function TZx(y) with the chosen normalization 
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FIG. 4. (Color online) Sketch of the numerical procedure for 
the 1 + 1 continuous DP in the finite box [0,tm] x [—ym,yrn] 
with the initial condition at to = 0.1 taken purely thermal 
V{to,y) =Vth.{to,y) (red curve). There are three levels of 
discretization grids by decreasing scale: the grid for the 
random potential ^f"^"* — 1 and ^f"^"* — 2 (the grey box on 
the left side); the linear grids for the recording of data 
At"" = tm/80 and Ay"" = 2j/„/100, and the logarithmic grid 



log 



to • (imAo)'^*° with i = 0.. 



the microscopic grid 



t 

for the numerical integration adapted in both {t, 3/)-directions 
by Mathematica. At a given configuration of disorder V{t,y) 
as illustrated in Fig. 22, the PDF 'Pv{t,y) is computed at 
increasing 'time' following (26), with the boundary condition 
riv{t,±ym) = and the ad hoc normalization (87), as illus- 
trated at to <ti < t2 (blue curves). 



/jjj dxTZ^{x) = 1. Firstly TZ^{y) is taken to be a Gaussian 
function, whose single feature is given by its variance 2^^: 



7^' 



Gauss 



(y) = 



4< 



Secondly we add phenomenologically two negative oscil- 
lations of the observed correlator, by using a cardinal 
sinus of period ^/tt with a Gaussian envelope function: 



SincG 



(y) 



sin Try/S, 



ny Erf(7r) 



(89) 



And thirdly we consider the same function as the exact 
two-point correlator of the microscopic random potential 



-^CubicS(y) 



^CubicS(j^) defined in (I2)-(I3). {|t. A} are 



reliable quantities if they do not depend on the choice 
of the fitting function, except for a numerical constant 
between Te*^^"''", 7^S'"=^ and n'^^^'''^ . I n Appendix I, 
we have compared the fit of the correlator V{t, y)V{0, 0) 
with respect to 7^Gauss^ ^SincG ^CubicS^ g^^^ deter- 
mined the numerical constants for passing from one to 
the others, as a consistency check of this procedure on 
this well-controlled correlator. 



B. 'Time'-evolution of the effective disorder 
correlator R{t, y) at fixed temperature 

In this section we study in detail the 'time'-evolution of 
the correlator R{t, y) measured numerically for three dif- 
ferent temperatures T G {0.35, 1, 6}, respectively at low-, 
intermediate- and high-temperature - cf. Appendix H for 
the complete set of corresponding numerical parameters. 

In Fig. 5, we can follow the evolution of the correla- 
tor, starting by construction from the thermal condition 
R{to, y) = at initial 'time' to = 0.1. At small 'times' the 
central peak first increases but quickly saturates, and all 
the curves start accumulating in the vicinity of y = {left 
side). Assuming thus that at sufficiently large 'times' the 
correlator has reached its presumably stationary form at 
small y, the saturation correlator i?sat(?/) can be obtained 
by averaging the correlator over t > tj^in (superimposed 
black curve in the right side). The corresponding satu- 
ration correlator ^sat(y) can then be fitted with respect 
to the three fitting functions (center), the corresponding 

values for the fitting parameters |fsat, i^sat| being listed 
in Tab. I. Having checked that those results are stable 
for different values imin > 10, we have chosen arbitrar- 
ily ^min = 25 for all temperatures in order to be safely 
above the saturation 'time' tgat which clearly exists but 
might be shifted by the artificial thermal condition at 
initial 'time' tg = 0.1. The negative excursions of R{t,y) 
at large y are treated as noise for the averaged Rgatiy), 
which as an artifact slightly displaces the asymptotic be- 
havior below the abscissa axis. 

So we can distinguish two regimes for the correla- 
tor R{t,y): a short-'time' evolution a priori marked by 
the artificial initial condition, and above a certain 'time' 
^sat ~ 10 a saturation of the central peak to a given sta- 
ble function Rsat{y) ~ -Dsat • ^(y) which can then be 
measured numerically provided that we neglect the neg- 
ative excursions at large y as proposed in our DP toy 
model (44). 

Going one step further and testing different fitting 
functions for TZ{y), we observe graphically that at high- 
T the exact microscopic disorder correlator TZ''^^^^'^^ (y) or 
alternatively its phenomenological counterpart T^SmcG 
both correctly encompass the features of the whole peak, 
including its maximum and its negative anti-correlations, 
but not the negative asymptote which is an expected arti- 
fact of the averaging procedure. This behavior could thus 
be consistent with the decomposition (42)- (43) at high-T 
and at 'times' t > isat, and the infinite-'time' limit (39) 
where the microscopic disorder correlator is recovered. 
At low-T however TZ'^'^^^'^ (y) appears to be more suited 
to capture the central peak and its maximum, although 
completely skipping the anti-correlations which are obvi- 
ously inherited from the microscopic disorder correlator 
but progressively altered when the temperature is low- 
ered. 

Assuming that the validity of the decomposition 
R{t, y) Dt ■ TZt (y) could be extended to smaller 
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r = 0.35 



2.C 




R{t,y) -Rsat(y) 





FIG. 5. (Color online) Effective disorder correlator R{t, y) 
measured numerically at fixed temperature T £ {0.35,1,6}. 
The different 'times' t£ [0.1,40] are separated in two sub- 
sets [0, tmin] {left) and [tmin,40] {right) with tniin = 25, with 
a timestep At = 0.5. Left side {y £ [-20,20]): Initial de- 
velopment of the central peak, saturation and accumula- 
tion of the curves in the vicinity of y = 0. Right side 
(j/ € [—40,40]): Large 'times' correlators and their average 
over t > tmin with tmin = 25 (superimposed black curve). 
Center: Saturation correlator Rsa.t{y) (black dots) with its 
Gaussian (88) (red) and 'SincG' (89) (blue) fitting functions, 
'CubicS' collapsing exactly on 'SincG' - cf. Tab. I for the 
explicit values of the fitting parameters. 



{t > tmin = 25) 


T = 0.35 


T = 1 


T = 6 


T-)Gauss 


4.6 ±0.2 


4.3 ±0.2 


1.13 ±0.04 


T-)SincG 


4.4 ±0.2 


3.9 ±0.2 


0.99 ±0.03 


p^CubicS 
-^sat 


4.37 ±0.07 


3.89 ±0.05 


0.98 ±0.06 


T-)SincG / r~)Gauss 
J-^sat 1 -^sat 


0.95 


0.91 


0.87 


^Gauss 
Ssat 


0.61 ±0.04 


0.70 ±0.04 


0.92 ±0.04 


^SincG 
Ssat 


2.25 ±0.07 


2.42 ±0.05 


2.99 ±0.06 


P'CubicS 
Ssat 


2.25 ±0.07 


2.40 ±0.05 


2.97 ±0.06 


P'SincG /^Gauss 
Ssat / Ssat 


3.72 


3.45 


3.25 



TABLE I. Fitting parameters {^sat, -Dsat} at T G {0.35, 1, 6}, 
computed from the saturation correlators Raa.t{y) plotted in 
Fig. 5 for tmin = 25. 



'times', we have followed the evolution of the two fitting 
parameters {^t,£)(}, as plotted in Fig. 6-7 for the same 
fixed temperatures T S {0.35,1,6} as in Fig. 5. Even 
though this assumption a priori breaks down for t < tsat, 
the three fits yield consistent values both for Dt and ^t, 
with reasonable uncertainties even at short-'times'. On 
one hand we recover for Dt the two regimes observed for 
the correlator R{t,y), i.e. an increase of this amplitude 
at short 'times' and a saturation beyond tsat as indicated 
by the straight lines _Doo in Fig- 6 (obtained by averaging 
Dt over t > tmin = 25) . On the other hand slightly de- 
creases at short 'times' and its saturation seems to appear 
much sooner, especially at high-T. As a self-consistency 
check, we can notice that the average of the fitting pa- 
rameters |^oo,i^oo| have the same values as the fitting 

parameters of the averaged correlator |^sat,-Dsat| (cf. 

Tab. I); their estimated errors are actually surprisingly 
close even though they have distinct origins, stemming 

respectively from the variance of Dt | at i > taun and 

from the uncertainty over the fit of Rsntiu). 

If the function TZ^^{y) did coincide exactly with 

the microscopic disorder correlator i?^"'^"^^(y) at 

we would expect that 

« 3.6 as dis- 



small y and large 'times' 



0.80 and 



nGauas £)Gauss — cGauss f 

cussed in Appendix I. Although |^oo,£'oo| clearly ex- 
hibit a temperature dependence which will be discussed 
later on, these ratios are numerically found to be of the 
expected order of magnitude and reasonably constant for 
t > isat at fixed T (cf. Tab. I). Actually the numerical 
ratios extracted from the measured R{t, y) versus the lin- 
earized i?''"(<,j/) of Appendix E (for i?C"bicS(y)) ^^q^ 
qualitatively the same trend in temperature, and are even 
quantitatively the same at high-T. This last point sup- 
ports furthermore the idea that at high-T we recover the 



microscopic disorder correlator R{t,y) « ^ • TZ^ '^°{y) 

but the average over t > tmin introduces an overall neg- 
ative shift of the asymptote which eventually yields the 
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T = 0.35 J 



rjGauss 

rjSincG 
~ °° 
rjCubicS 

rjSincG / r) Gauss 
10 



4.6 ±0.3 

4.4 ±0.3 
4.4 ±0.3 
0.95 



T = 0.35 



Pi 



k 



cGauss 

tSincG 
Soo 

cCubicS 
t-SincG /cGauss 



0.61 ±0.02 
2.25 ±0.05 
2.25 ±0.05 
3.72 



fa 



'fTBmHt| miT 5iTTTlWnT(tgi. ni TrnHTm r pihiTtinTniTTinr |ft» £ 



D, 



p8 




^ Gauss 
rSSincG 
riCubicS 



4.3 ± 0.2 
3.9 ±0.2 
3.9 ± 0.2 
0.91 



T = 1 



T= 1 



W[Blii i iffl|ffi%aHnii,,ii#iiBiiitf'giiii'%,i 



cGauss 
Soo 

tSincG 
Soo 

cCubicS 

Soo 

(^SincG /tGauss 
Soo / Soo 



0.70 ±0.02 
2.42 ±0.05 
2.41 ±0.05 
3.45 



0.6 
0.4 
02 



rjGauss 

-^oo 

^CubicS 
riSincG / rjGauss 



= 1.12 ±0.05 
= 0.98 ±0.05 
= 0.98 ±0.05 
0.87 



r = 6 



FIG. 6. (Color online) Fitting parameter Dt at fixed tem- 
perature T £ {0.35,1,6} with respect to 'Gauss' (red) and 
'SincG' (blue), 'CubicS' collapsing exactly on 'SincG'. The 
straight lines indicate the respective average value Daa for 
t > imin = 25, whose standard deviation is given explicitly for 
each fit. 



2.5 J 



cGauss 
Soo 


= 0.92 ±0.02 


cSincG 

Soo 


= 2.99 ±0.09 


cCubicS 
Soo 


= 2.97 ±0.09 


SincG /cGauss 

oo / Soo 


3.25 



FIG. 7. (Color online) Fitting parameter at fixed tem- 
perature T £ {0.35, 1, 6}, computed conjointly to Dt of Fig. 6 
and following the same convention regarding the legend. 



fixed ratios sa 0.89 and w 3.27. 

Beyond tiiese asymptotic values of the fitting parame- 
ters, their short-'times' behavior can have either a deep 
physical meaning or be an artifact due to the inadequacy 
of the fitting function TZ-iy) for the measured correla- 
tor R{t,y). The short-'times' increase of Dt will be dis- 



cussed from an analytical point of view in Sec. VB, but 
for now we compare in Fig. 8 the evolution of the peak 
^•(^,0) deduced from the fits (combining 7e^^^(0), Df*- 

and ^f*) to the values measured numerically as a direct 
consistency check. In Fig. 8 (top) we see that at low-T 
the peak is correctly captured by 'Gauss', and increas- 
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rRit,0) 



T = 0.35 




J?(t,0) 




Ai?,,t(0)/i?,,t(0) 



FIG. 8. (Color online) Top: 'Time'-evolution 

of the maximum of the correlator peak R{t, 0) at 
fixed temperature T £ {0.35, 1, 6}, measured numerically 
(black dots) or predicted by the fitting functions via 
i?^'(t,0) = 7^|ll(0) ■ iPfVlf (red for 'Gauss' (88), blue for 
'SincG' (89), cyan for 'CubicS' (I2)-(I3)). The curves 
for 'SincG' and 'CubicS' systematically coincide. Bottom 
left: Zoom on the high-T case, where the 'SincG' and 'Cu- 
bicS' predict more accurately the peak at high-T, whereas 
at low-T 'Gauss' seems more suited. Bottom right: De- 
crease in the relative differences between the 'Gauss' and 
'SincG' parameters with increasing temperature, averaged on 

A ^SincG ^Gauss 

the large-'time' window t £ [25,40]: - °° = °° ~sincG° > 



A52C 

and 



(cf. end of the Appendix. I) 



Afloat (0) 



"(o)-R°r°'(o) 



'(0) 



ing the temperature the relative difference between the 
fits decreases anyway; however, zooming on the high-T 
{bottom left) we see that 'CubicS' and 'SincG' catches 
more precisely the peak. Comparing finally the rela- 
tive differences between the fits as a function of T {bot- 
tom right), two temperature regimes can be postulated 
regarding the normalized function TZ{y)'. at high-T we 
have TZ{y) « 7?.'-""*'"^^ (y) with the relative behavior be- 
tween the fits as expected from Appendix 1 (geometrical 
collapse of ^^"^G/I^auss ^ ^SincG/^Gauss^ overcstimation 
of the amplitude and thus of the peak by 'Gauss'); at 
low-T the modification of the function Tl{y) is essentially 
pushed into the increasing discrepancy A^oo/'Coo- 

As a last consistency check before discussing the ac- 
tual temperature-dependence of the parameters ^oo and 
Doo, we have measured the evolution of the integral 
2 J^"^ dyR{t, y) that should be exactly zero if ym oo. 
On our restricted numerical window however this integral 



fiuctuates indeed around zero but very noisily, and once 
normalized by Doo it remains bounded by 1 consistently 
with respect to the decomposition (42)- (43). 



C. Temperature-dependence of the asymptotic 
effective disorder RsRt{y) 

In the previous section we have illustrated the 'time'- 
evolution of the correlator R{t, y) for three typical tem- 
peratures, and discussed its short-'times' versus satura- 
tion and low- versus high-T regimes regarding the shape 
TZ{y) of its central peak. Now we aim to characterize 
the temperature-dependence of the amplitude Doo and 
typical spread ^oo of the correlator, and point out the 
crossover between the two temperature regimes over the 
range T g [0.35, 6] - cf. Appendix H for the complete set 
of corresponding parameters. 

In Fig. 9 we have juxtaposed all the temperature- 
dependent Dt, £,t and Rsat{y)- We observe graphically 
that the amplitude Dt decreases strongly with T and the 
curves tend to accumulate at lower T, whereas for the 
typical spread £,t, it slightly increases (up to 30% never- 
theless). These behaviors result in an overall damping of 
the effective disorder correlator ^sat {y) when the thermal 
fluctuations are developing. 

The measured amplitude Doo as a function of T 
is reported in Fig. 10 for both 'Gauss' and 'SincG' 
(which collapses with 'CubicS') fits. Qualitatively 
it decreases in 1/T at high-T and saturates at low- 
T, accordingly to the scaling arguments and GVM 
predictions of the full model and the DP toymodel 
(cf. Sec. IIIB-IIIC). A parametrization of this tem- 
perature crossover has been defined in the relation 
Doo(T, ^) = /(T, ^) (45) with the interpolating param- 
eter /(T, ^). Quantitatively its ^ = limit requires 
at high-T that / ^ 1 without any additional numeri- 
cal prefactor, but extracting the strength of disorder D 
from the 1/T behavior of Doo actually yields a system- 
atic underestimation with respect to the microscopic dis- 
order: T»Gauss ^ g j2 ± 0.08, dS'"<=G = 5.92 ± 0.05 and 
^CubicS ^ 5 gg _|_ Q_Q5 (obtained on the three larger avail- 
able temperatures), whereas D — 8 (86). We attribute 
again this discrepancy to the negative excursions of 
R{t, y) at large y which bias all the fits and preclude 
a quantitative test of /(T, ^). 

The crossover temperature cannot be sharply de- 
termined from Doo(T, ^), nevertheless the definition 
Tc{0 = (CcD)i/3 obtained by GVM (46) and by scahng 
arguments (56) predicts with ^ = 2, c = 1 and D = 8 that 
Tc « 2.5. Without any corrective numerical constant this 
value is actually compatible with the crossover of the am- 
plitude in Fig. 10 and also of the function Tl{y) in Fig. 8 
{bottom left). 

Finally the typical spread ^oo obtained in parallel to 
Doo is reported in Fig. 11. The collapse of the different 
fits for ^oo has already been discussed in the previous 
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FIG. 10. (Color online) Temperature dependence of Doo with 
respect to 'Gauss' (red) and 'SincG' (blue). The straight line 
cD/T indicates the expected behavior in the high-T regime 
(cD = 8 for our data). Inset: Corresponding mean and stan- 
dard deviation for the average of Dt over t > tmin = 25. 



FIG. 9. (Color online) Temperature-dependent Dt and 
for 'SincG' {top) and Rsatiy) with tmin = 25 (bottom), for 
T G [0.35, 6] as listed in Appendix H. For increasing tempera- 
tures (blue to red), Dt decreases whereas slightly increases 
|-^grid _ 2 fQj. reference), resulting in an overall damping of the 
correlator Rsat{y)- The small negative excursion at low-T for 
is an example of the numerical sensitivity of the noidinear 
fitting procedure. 



subsection; it is compatible with the scenario of a high-T 
correlator 7?.sat(2/) ~ T?.*""*^"^^ (y) . However it displays a 
temperature-dependence that a priori corrects to first- 
order the minimal assumption that « ^ in our DP 
toymodel; since the thermal fluctuations seem to increase 
the effective ^oo compared to the microscopic disorder 
correlation length ^siid = 2 (which consistently remains 
a lower bound in our measured '^^"'^"^^). We know from 
the linearized solution (40) that neglecting the KPZ non- 



linearity we recover asymptotically %j^TZ' 



cU-nCubicS 



(y) for the 



correlator, so any modiflcation of in this two-point 
correlator can only stem from the KPZ nonlinearity at 
high-T. A numerical artifact similar to the underestima- 
tion of Doo is not to be excluded, but no more conclusions 
can be drawn from our numerical results. 



D. Fluctuations of Fv{t,y) 
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FIG. 11. (Color online) Top: Temperature dependence of 
with respect to 'Gauss' (red) and 'SincG' (blue), normalized 
for 'Gauss' with ^SincG^^Gauss _ 3 g ^^^^ Appendix I). Bot- 
tom: Corresponding mean and standard deviation for the av- 
erage of over t > tmin — 25. 



In its initial definition (12-13), the disorder free-energy 
Fy(t,y) is defined up to a constant ctev'(i) depending 
of the chosen path-integral normalization of Wv(t,y) in 
(3). As far as statistical averages with quenched dis- 
order are concerned, such as (6) and (22), this con- 
stant is irrelevant and thus usually completely skipped. 



The Feynman-Kac equation (26) actually yields a uni- 
vocal definition of Fv{t,y) which satisfies the STS (17), 
and as a consequence the evolution of its y-independent 
mean value is given by (33) starting from the ini- 
tial condition Fv{0,y) = 0. In our numerical approach 
Fv{t, y) = cte(t) is in fact a tractable quantity which pro- 
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FIG. 12. (Color online) Disorder average of Fv{t,y) at 
fixed temperature T G {0.35, 1,6}. Top: Fv{t,y) as a func- 
tion of y, for increasing 'times' f G [0.1,40] with 'time'-steps 

At = 2.5 (top to bottom). Bottom: —2cFv{t,y) averaged 
over 1/ as a function 'time', the error-bars indicating the cor- 
responding standard-deviation; the result of the linear fit for 
t > tmin = 25 is indicated by the straight lines in black. 



vides an independent way to measure Rsntiy = 0) from 
its large-'times' linear evolution. This is remarkable and 
a direct consequence of the KPZ nonlinearity, since lin- 
earizing (26) trivially predicts F|j"(i,?/) = 0. 

In Fig. 12 we have plotted the disorder average of 
Fv{t^y) at our usual fixed temperature T e {0.35,1,6}, 
keeping first the spatial resolution in [t, y) {top) and then 
averaging over the ^/-direction {bottom). Fv{t,y) should 
be y-independent, but at short-'time' and low-T it dis- 
plays a slight curvature that we attribute to the artifi- 
cial thermal condition at to = 0-1; the resulting standard- 
deviation is strongly reduced at higher T where this ini- 
tial condition is more accurate with respect to thermal 
fluctuations. 

At large 'times' Fv{t,y) follows a robust linear be- 
havior which extends down to tsat ~ 10, as emphasized 
in Fig. 12 and generalized to all the temperature spec- 
trum in Fig. 13. Below tgat no clear powerlaw could be 
identified though the logarithmic scale in Fig. 13 makes 
explicit a superlinear short-'time' behavior, a priori con- 
ditioned by the initial thermal condition but see Sec. V B. 

We observe graphically in Fig. 14 that the slope is 
damped with increasing T, which corresponds to the de- 
crease of the maximum Rsa,t{y = 0) and physically to the 
damping of the effective disorder due to larger thermal 
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FIG. 13. (Color online) Temperature-dependent 

— 2cFv{t,y) averaged over j/ as a function of 'time', for 
T G [0.35, 6] as listed in Appendix H. For increasing temper- 
atures {blue to red) the slope decreases as expected by (33). 
Inset: Zoom on the short-'times' behavior in a logarithmic- 
scale, which shows qualitatively the existence of a saturation 
'time' tsat ~ 10 marking the beginning of the linear regime at 
large 'times' (the dotted lines guide the eye for a quadratic 
behavior ~ t^). 



fluctuations. A quantitative comparison of the different 
values of .Rsat(O) function of T is given in Fig. 14, 
including a direct measurement on i?sat(y), the predic- 
tions obtained by the flts, and the independent mea- 
surement via Fv{t,y). We find an excellent agreement 
between these quantities, with a slight crossover of the 
direct measurement from the 'Gauss' to the 'SincG' fits, 
and a systematic overestimation of —2cdtFv{t,y)', this 
last point is an artifact of the short-'times' curvature in 
y oiFv{t,y). 

After this successful consistency check, one could won- 
der what would be the quantity which is the more rele- 
vant physically: the maximum Rsa.t{y — 0) or the ampli- 
tude I?oo? The first can be measured independently from 
dyFv{t,y)dyFv{t,0) and Fv{t,y), whereas the second is 
obtained assuming a given fitting function TZ{y) and bi- 
ased by the negative excursions at large y (the 'wings' 
of C{t,y)). In the limit ^ = on one hand ^sat(O) di- 
verges, whereas Dao remains finite, is given by ^ and its 
zero-temperature limit is thus ill-defined. At finite ^ > 
on the other hand i?sat(0) remains finite, tends to zero 
at high-T and saturates at low-T as does its amplitude 
Doo ■ So we need a finite ^ > in order to define _Rsat (0) 
and tackle properly its different limits in temperature; 
Doo is always regularized and captures the physical con- 
tent of -Rsat(O)- being thus better suited to characterize 
a potential universal behavior. 
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FIG. 14. (Color online) Temperature dependence of the max- 
imum of Rs!i.t{y) over 'times' t G [25,40], measured by three 
different ways: first measured numerically directly .Rsat(O) 
(black crosses), secondly deduced by the 'Gauss' and 'SincG' 
fits as in Fig. 8 (respectively red and blue dots), and thirdly 
measured indirectly from the linear slope of —2cFv{t,y) in 
Fig. 12 (purple dots) and systematically slightly overesti- 
mated. Inset: The excellent linear behavior of Fv{t,y) yields 
vanishing error bars for the slope on the range t £ [25,40]. 



E. Temperature-dependence of the roughness 
function B{t,T) and its exponent (^{t,T) 

Up to now we have analyzed in detail the disorder 
free-energy fluctuations, specifically the two-point corre- 
lator R{t,y) = ^dyC{t,y) and the mean Fv{t,y). Keep- 
ing in mind the equivalence between the 'time' t for 
the DP endpoint and the lengthscale r on the ID in- 
terface (cf. Fig. 1), we focus now on their consequences 
on the geometrical fluctuations, and speciflcally on their 
variance as a function of 'time', namely the roughness 

As recalled in Sec. II C after the deflnition (11), the 
roughness is expected to display at least two asymptotic 
powerlaw regimes, crossing from a pure thermal behav- 
ior B{t) > i3th(i) = ^-t short 'times' to the random- 
manifold Basympt(i) ~ D^t'^l'^ at large 'times' - as 
hinted by scaling arguments (66)- (67) and thus also pre- 
dicted by GVM (50) . Above there is a single crossover 
Larkin 'time' t*(T) = (54) which diverges with T, 
and essentially prevents us from observing the completed 
crossover to the RM regime in our numerical approach in 
a reasonable computational time for T > 2. Below Tc 
however we expect from Ref. 34 the appearance of an in- 
termediate 'Larkin-modifled' roughness regime ending at 
the generalized Larkin 'time' Lc(T^^) (49) which marks 
the beginning of the RM regime . 

Though the GVM framework provides us a prediction 
^DP (t) for the full temperature-crossover of our DP toy- 
model (cf. Appendix A), a complete quantitative test of 
(A7-A10), (A16) and (A17) is hindered by the numeri- 



cal constants that are a priori present in the determi- 
nation oi D, ^, Tc, Lc and the amplitude of the asymp- 
totic roughness itself A{c, D,T,^), due to previously dis- 
cussed numerical artifacts and also to the GVM approx- 
imation. Nevertheless, we observe qualitatively these dif- 
ferent roughness regimes, as plotted at fixed tempera- 
ture T S {0.35, 1, 1.8} in Fig. 15, indicating both the to- 
tal roughness B{t) — {y{t)^) and its pure disorder com- 
ponent Bdis{t) = B{t) - ^ = {yi^Y) as stated by the 
STS (B12)-(B13). In all three graphs we can follow 
the crossover from the thermal asymptote (red) to the 
RM asymptote (green) which stems from the increasing 
Bdis{t) added to Bth{t). Note that the RM asymptote 
has been constructed consistently with (50) with a nu- 
merical correction fixed once and for all from the dataset 
T = 0.4 averaged over t G [25, 40], assumed to be already 
in the RM regime: 
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0.292 ± 0.008. 



(90) 

This RM asymptote 



is graphically consistent with all the available datasets 
in the range T e [0.35, 1.8], which plays in favour of a 
numerical corrective factor common to all temperatures 
and absorbing the discrepancy in i3asympt(^) stemming 
from Doo and the GVM. The low-T regime is illustrated 
by r = 0.35 where the intermediate 'Larkin-modified' 
regime is clearly present, whereas it has disappeared as 
such already at T = 1. 

Larger 'times' have been explored at T = 1.8 which is 
believed to be close to (cf. Sec.IVC), and is consistent 
with lower temperatures except for a sudden increase at 
t > 400 which can clearly be attributed to the finite size 
of the box y £ [— ?/m,?/m] ^ the condition B[t) < y^ 
translates into B{t) < 160^ = 2.5 • 10* for this dataset 
and yields an adequate upper bound in Fig. 15 (bottom). 
The corresponding upper bound for most of the datasets 
is B{t) < 40^ = 160 thus the range of roughness displayed 
in all the other graphs. 

Gathering in Fig. 16 the roughness B{t) and i?dis(0 
over the range Te [0.35,6], we observe as expected 
with increasing temperature that the disorder rough- 
ness B^isit) is progressively damped by thermal fluctua- 
tions and thus that the intermediate 'Larkin-modified' 
regime shrinks with increasing temperature. The be- 
ginning of the RM regime is consequently pushed to 
larger 'times', as it could also be deduced by the con- 
dition that Bdis{t) becomes comparable to Bth{t) close 
to the Larkin 'time' t — Lc- From the point of view of 
the ID interface, the temperature-dependence of i?dis(0 
can be physically understood with the following picture: 
at small lengthscales the thermal fluctuations make the 
interface rougher within the local valleys of the disor- 
dered free-energy landscape and B{t) > Bth{t) increases 
with T as expected; at large lengthscales on the contrary 
thermal fluctuations allow the interface to explore more 
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FIG. 15. (Color online) Roughness B{t) = {y{ty^) and 
its disorder component Bdis{t) = B{t) — = (vi't)'^) mea- 
sured numerically at fixed temperature T £ {0.35, 1, 1.8}. 
Main: Focus on the crossover of B{t) (thick) from the pure 
thermal behavior Bth{t) = ^ (dotted red) at short 'times' 

to the asymptotic Basympt(t) ~ Doo{T)^^^f^^''^ (dotted green) 
at asymptotically large 'times'; Bdis{t) is the thinner lower 
curve. Inset: Focus on the short-'times' crossover of -Bdis(i) 
(thick), in parallel to B{t) (thin). 



effectively the free-energy landscape by overcoming some 
free-energy barriers in order to minimize its elastic en- 
ergy and thus have i?asympt(0 decreasing with T. From 
the point of view of the DP, these behaviors are encoded 
in the evolution of the translation-invariant distribution 
V [F,t\ (17), the integration of the microscopic disor- 
der V{t' , y') explored by the elastic DP over 'times' [0, t] 
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FIG. 16. (Color online) Temperature-dependent rough- 
ness B{t) = {y{tY) and its disorder component Bd^it) (in 
logarithmic scale, respectively top and bottom curves), for 
T G [0.35, 6] as listed in Appendix H (blue to red for increas- 
ing temperatures). Bdiaif) for T = 6 is indicated in dotted 
red. 



defining the effective disorder Fy{t,y) and consequently 
the DP-endpoint fluctuations (22). The saturation at 
low-T of the amplitude Doo in Fig. 10 and of the peak 
R{t, 0) in Fig. 14 naturally translates into a saturation of 
the asymptotic amplitude of the roughness (90). 

In order to focus on the disorder contribution, we 
have computed numerically directly Fv(t, y), but regard- 
ing the roughness we had first to construct the total 
free-energy F = Fth + F, then compute the total rough- 
ness B{t) and eventually deduce its disorder compo- 
nent i?dis = B — Bth- This quantity is thus more sub- 
ject to noise at higher T (see i.e. T = 6 in Fig. 16), 
but especially at low-T it clearly displays two regimes in 
'time'. To address the question of a possible powerlaw at 

short-'times', the logarithmic slope C(dis)(i) — \ ^^aiogt^ 
is plotted in Fig. 17. While C(i) crosses over as ex- 
pected from QtYi = 1 to Cbm = | but excludes the def- 
inition of an intermediate-'times' powerlaw for B{t), 
Cdis(i) on the contrary displays a plateau at low-T which 
disappears already at T = 1, a tendency confirmed at 
T = 1.8. Accordingly to the GVM prediction (A8) 

Bdis{t) should start in ^ H in which case the value 
of the plateau would have 1 as a lower band and only a 
'time'-dependence of Dt and/or could correct Cdis to 
match the GVM prediction, which is anyway not to be 
trusted a "priori at 'times' shorter than Lc (cf. Sec. IIIB). 
The value of this plateau exhibits moreover a slight 
temperature-dependence, which cannot be accounted for 
by our DP toymodel but might also simply be an arti- 
fact of the thermal condition imposed at — 0.1. Note 
finally that all these effects on i?dis(i) take place before 
^sat ~ 10 as estimated on the saturation of the free-energy 
correlator R{t, y) in the previous sections. 

As a last remark in this numerical study of the 1+1 DP 
fiuctuations with a finite disorder correlation length, we 
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FIG. 17. (Color online) Logarithmic slope of the roughness, 
("(t) and (^dis(i) obtained respectively from B{t) and Bdis(t) 
of Fig. 16, for T £ [0.35, 1] (blue to red for increasing temper- 
atures) and T — 1.8 (in dashed dotted black). The thermal 
exponent C,th ~ \ and the RM exponent ("rm = | are indi- 
cated by the dashed dotted lines, respectively red and green. 
At larger 'times' those logarithmic slopes are very noisy but 
still oscillate around the expected value C,i{.m = I. 



check the assertion that the PDF V{t, y) of the DP end- 
point is 'fairly Gaussian' (10), already checked at f = 
in the RM asymptotic regime'*'^. We have computed in- 
dependently the higher moments {\y{t)\^) (2 < fc < 6) of 
the PDF and renormalized them by the Gaussian 'time'- 

independent coefficients \/2 [F {^^) /-v/tt] ^^'^ with F the 
Euler function. If V{t,y) is a pure Gaussian function, 
these renormalized moments should perfectly collapse on 
each other at all 'times', as it is the case in at short 'times' 
where V ~ 'Pt\i and remains reasonably so by naked-eye 
in the intermediate 'Larkin-modified' regime. At 'times' 
above the Larkin 'time' we observe only a very light fan- 
ning of the renormalized moments, showing the consis- 
tency of our observations at f > with respect to the 
previous numerical studies at ^ = 0. 



V. DISCUSSION 

A. Synthetical outlook 

We addressed throughout this paper the consequences 
of a finite correlation length ^ > of the microscopic dis- 
order V{t,y) explored by a 1+1 DP, or alternatively of 
a ID-interface finite width which is always present in 
experimental systems. On one hand, several analytical 
arguments yielding exact results at ^ = break down 
as such, questioning their generalization to ^ > 0. On 
the other hand, despite a lack of exact analytical ex- 
pressions the finiteness of this quantity allows to control 
the scalings and the low-temperature limit of the model 
(see Sec. IIIC-IIID), avoiding the pathological and un- 



physical divergences that appear at ^ = 0, in particular 
conjointly to the limit T — > 0. 

In order to tackle the case at ^ > 0, the 1+1 DP formu- 
lation allows to follow effective quantities at fixed length- 
scale or growing 'time', in an approach thus conceptually 
similar to the FRG which focuses on the flow and fixed 
points of the disorder correlator (denoted A or /j32,33,64-j^ 
Considering the free-energy at fixed disorder (averaging 
over the thermal fluctuations but one step before the dis- 
order average), it is thus possible to disconnect theoret- 
ically and numerically the two statistical averages, and 
even to focus on the pure disorder contributions thanks to 
the STS and the Feynman-Kac equations for Fy and its 
derivative r/y (see Sec. II D). Those two quantities are not 
directly accessible experimentally (except for liquid crys- 
tals, as discussed later in Sec. VI) and are a priori more 
complex to handle since they encode more information 
than direct observables such as the geometrical fluctua- 
tions and the roughness. However they actually display a 
simpler phenomenology by disconnecting the thermal and 
disorder effects and the different lengthscales, whereas 
the roughness B{t) intricates all of them. 

Although V \F\ and V [rj\ are not Gaussian at ^ > 
not even in the inflnite-'time' limit, their main fea- 
tures are encoded in the two-point correlators C{t^ y) 
and R{t,y), i.e. the scalings of C and R dominate the 
higher moments of the PDFs, similarly to the case f = 
(see Sec. Ill A). This supports consequently the construc- 
tion of the toymodel of Sec. IIIB, which relies on the as- 
sumption that the PDFs can be approximated as Gaus- 
sian ones described by a given set of two-point correla- 
tors. The numerical results of Sec. IV are indeed quali- 
tatively in agreement with the GVM predictions derived 
from this DP toymodel (see Appendix A) , in the range of 
numerical parameters that we have considered (see Ap- 
pendix H) . The study of the two-point correlator R pro- 
vides thus a vantage point on the DP properties, flrst 
at asymptotically large 'times' (keeping in mind that the 
inflnite-'time' limit simplifles the analytical treatment of 
the RM regime, i.e. via a Fokker-Planck approach as 
in Appendix F) and secondly at flnite 'time' with the 
connection to the short-'times' regime. 

In order to characterize the asymptotic large-'times' 
behavior, which can potentially display universality, a 
central quantity is the saturation amplitude of the effec- 
tive disorder i.e. Dao{T, ^). At flxed ^ > it is equivalent 
to the maximum of the saturation correlator i?sat (y = 0) 
which is equal to -Doo/Coo • "T^i^ily = 0). However 
for a (5-correlated microscopic disorder R^=Q{y) = S{y), 

Rsatiy = 0) diverges whereas the amplitude I?oo/^oo re- 
mains well-deflned. It is remarkable to notice that from 
the whole saturation correlator, the quantity Doo is the 
only feature that eventually plays a role in the asymp- 
totic roughness in GVM or scaling arguments e.g. in (50), 
the speciflcity of the (normalized) RB disorder correla- 
tor R^{y) (1) and thus of the function TZ{y) being then 
gathered into a numerical constant. Doo(T,^) appears to 
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be the relevant quantity for a universal description of the 
crossover between low- and high-T asymptotic DP fluc- 
tuations, from an analytical point of view and in agree- 
ment with the numerical results in Sec. IV. Going one 
step further, the crossover from its ^ = (or high-T) 
limit is better described by the interpolating parame- 
ter f{T,£,) = l)oo(r,0/ (^) first introduced in (45) and 
expected to rescale also the characteristic scales such as 
Lc{T,C) according to the relations (58)-(61) obtained by 
pure scaling arguments. The GVM framework yields the 
two predictions (46) and (70) for /(T, ^) derived from the 
value of the fuU-RSB cutoff (cf . Appendix A) , and a third 
analytical prediction will be presented in the next sub- 
section VB; all of them predict a monotonous crossover 
connecting the limits T — and ^ = (62)-(64) with a 
polynomial equation on f{T,(,), which is qualitatively 
consistent with the numerical results (see Fig. 10) but 
suffers quantitatively from the variational approximation 
and several corrective factors as discussed throughout 
Sec. IV. Nevertheless an important and successful test 
for all this discussion is the independent but consistent 
numerical determinations of i?sat(2/ = 0) iii Fig- 14, also 
with respect to dtFv{t,y), and the adequacy of the de- 
duced asymptotes i?asyinpt(^) in Fig. 15, having system- 
atically tested and discussed the small and controlled 
discrepancies of the different fitting functions ('Gauss', 
'SincG' and 'CubicS' defined in Sec. IV A). 

Note that although the two communities of physi- 
cists and mathematicians work with two different con- 
ventions, respectively at fixed elastic constant c (the 
choice c — 1 essentially fixing the units of energy) ver- 
sus ai c — T (as discussed in Appendix G) , all the above 
discussion remains valid in both conventions, although 
the choice c = T leads to other limits in temperature. 
The two opposite limits at high-T (or f « 0) Doa ~ 
and at low-T (or ^ > and below T^iO = {£,cD)^/^) 
_Doo ~ — (c^D^^^^)^/-^ translate with the conven- 
tion c = T into Doc^ D and Dr^ « T2/3^2/3_e-i/3 
spectively above and below — \/£T) (deduced self- 
consistently from T^^ = {icDf/^ = {^T^Dy/^). In the 
course of the study of the rescaling of the correlator 
C{t,y) with respect to the roughness B{t) in Ref. 46, it 
has been noticed that in the regime \y\ < \/ B(t) we have 
numerically as expected a linear behavior C(t,y) cx \y\ 
but with a by-product prefactor that corresponds pre- 
cisely to our Doo ■ Taking as a criterion the collapse of the 
curves C{t,y) at different temperatures on an arbitrary 
chosen curve, the temperature dependence of this pref- 
actor is consistent with all our analysis on the origin and 
interpretation of Dao{T,^) (see the insets in Fig. 10 and 
Fig. 14 of Ref. 46, which illustrate respectively the con- 
ventions c = T versus independently fixed c = 1 and T) . 

As for the finite-'time' behavior, especially at short- 
'times' it is sensitive to the specific microscopic disor- 
der correlator, thus compromising a possible universal- 
ity. The 'time'-evolution displays essentially two regimes 
on the fluctuations of the disorder free-energy, separated 



by ^sat (whose order of magnitude has been crudely es- 
timated in our numerical study to tgat ^ 10): starting 
from the initial condition R{0, y) = 0, the central peak 
of the correlator R{t, y) develops itself keeping the inte- 
gral J^dy R{t,y) —0 constant, till it reaches the satu- 
ration shape T^^iv) compensated by negative bumps 
according to the generic decomposition (42)-j^43) as il- 
lustrated in Fig. 3. Since dyC(t,y) = 2R(i,y), after 
the double integration (21) the correlator C{t,y) starts 
from the initial condition ^(0, y) = and at fixed 'time' 
above isat it is rounded at \y\ < ^, increases then linearly 
C{t,y) « Doo\y\ at ^ < \y\y/B(t) and is constant for any 
larger \y\ (see again Fig. 3). On one hand, the position 
it of these 'wings' of C{t,y) or equivalently of the neg- 
ative bumps of R{t, y) is discussed at length in Ref. 46 
and identified to correspond physically to the typical po- 
sition of the DP endpoint, it ~ ^/B{t), in the different 
roughness regimes and even below the Larkin length Lc- 
On the other hand, what happens below tgat cannot be 
understood without taking into account the whole mi- 
croscopic disorder correlator R^{y); whose feedback via 
the KPZ nonlinearity at small \y\ modifies the ampli- 
tude Dao{T,£^) and the shape TZ^{y) above the satura- 
tion, especially in the low-T regime and in any case with 
^sat < Lc- Neglecting the KPZ nonlinearity yields the 
prediction (40)-(41) which mixes different limits: the 
^ = (high-T) amplitude I?oo = the same correla- 
tor as the microscopic disorder TZ^{y) — R^{y), and the 
'wings' rescaled with respect to the pure thermal rough- 
ness i3th(0 lengthscales (this diffusive behavior be- 
ing for sure an artifact of the linearization) . In the low- 
T regime, we believe that by generating relevant non- 
Gaussian correlations such as Rg, and C3 below tgat, the 
KPZ nonlinearity introduces an effective kernel for R that 
modifies simultaneously Doo and TZ^{y), with in partic- 
ular the saturation below Tc of the amplitude Doc, ~ ^ 
as predicted by scaling arguments in Sec. IIIC. 

The phenomenology of the DP fluctuations is simpler 
from the point of view of the disorder free-energy Fy, 
via its two-point correlators R and C which display two 
'time '-regimes separated by tsat- From the competition 
between the typical Fy and the thermal F^h the resulting 
roughness -Bdis(0 also displays two regimes as observed 
numerically in Fig. 16 and 17. However, when recom- 
bined with the pure thermal effect the roughness B{t) 
displays two or three 'time '-regimes respectively at high- 
T (thermal and RM regimes) and low-T (with an addi- 
tional intermediate 'Larkin- modified' regime), with Lc at 
the beginning of the RM regime, as predicted by GVM 
and observed numerically. The disorder free-energy is an 
effective quantity which encodes the microscopic disorder 
explored by the polymer, there is consequently a feed- 
back between the geometrical fluctuations V{t,y) and 
the free-energy correlations C{t,y): firstly the existence 
of 'wings' in C(t,y) are imposed physically by the finite 
variance of the PDF V{t,y) (the polymer does not ex- 
plore often regions |y| > ^JB{t) so these regions do not 
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contribute much to the correlator C{t,y)); secondly the 
PDF V{t^y) is deduced from the competition between 
Fth and Fv - C'^/^ the maximum of the typical Fy be- 
ing precisely fixed by the 'wings' of C] thirdly the 'wings' 
of C or the bumps in R can be skipped for a GVM compu- 
tation of the roughness, providing a self-consistent justi- 
fication of our DP toymodel. Beyond the scaling in 'time' 
of those fiuctuations, which we plainly understand phys- 
ically now, their temperature dependence at all 'times' 
is finally determined by the integrated disorder up to 
isat ^ Lc, where the KPZ non-linearity plays a crucial 
role below Tc{i). 



B. Effective evolution in 'time' and temperature of 
the amplitude Dt 

Having this global picture in mind, we can now gather 
all the physical intuition we have obtained and construct 
the following analytical argument in order to obtain an 
evolution equation for the amplitude Dt. 

The evolution of the correlator R{t,y), given by the 
'flow' equation (34), cannot be solved directly since 
it brings into play the three-point correlation function 
R3{t, y), a hallmark of the KPZ non-linearity. To extract 
an exact information from this flow one should in prin- 
ciple solve the full hierarchy of equations connecting the 
whole set of ?i-point correlation functions, a task which 
seems however out of reach. As we will detail there- 
after, the restriction of the flow to the vicinity of y = 
leads in fact to an (approximate) closed equation on the 
height R{t, 0) of the two-point correlator. It will allow to 
pinpoint the role of the non-linearity in the temperature- 
dependence of the asymptotic Doo{T,£^) and its interpo- 
lating parameter /(^,T) defined by (45), and give more 

insight into the short-'time' behavior of Dt and Fv{t,y) 
(with respect to (33)). 



1. Rescalmgs of R, R, R^ and -qy 
From (34), the flow of R{t, y) in y = reads 

dtR(t, 0) = -R"(t, 0) - -RUt, 0) - -R(t, 0) - DR'I(0) 
c c t ^ 

(91) 

(throughout this section we denote for short the deriva- 
tive with respect to j/ by a prime). Although this equa- 
tion is exact, it cannot be solved directly since the three- 
point correlator i?3 is not known. To go further and try 
to find out what relations between the physical param- 
eters it might nevertheless imply, one has to surmise a 
(minimal) scaling form of the different correlators and 
their first derivatives. 

Let's first consider the known scaling of the micro- 



scopic disorder correlator R^{y): 
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-R^ (0)1^=1 and ci = 

merical constants, independent of ^ and refiecting the 
specific geometry of the correlator around the origin. 
For instance when the correlator is a Gaussian function 
R^{y) ^nf '''''''' (y) defined by (88), one has co = ^ 

and ci = |. 

By analogy with (92) and supported by the numerical 
test of our DP toymodel in Sec. IV B, we now assume 
that the correlator R{t, y) scales around y « as 



(93) 



Here, C2 and C3 are numerical constants independent of 
the parameters (c, D,T,£^,t): 



C2 = i?(t, 0) 



Dt = l 



C2C3 



-R"{t,0) 



(94) 



while is the height of the central peak R{t,0) as- 

sumed to capture all the dependence in the parameters. 
C2 is actually defined so that the ^ ^ limit (38) is 
recovered: 

D^ (T, e) = lim Dt^ [ dy i?(oo, y) > (95) 
cD 

limZ?^(r,e) = — (96) 

which is known to hold exactly, without any additional 
numerical constant. The constant C2 depends on global 
properties of the infinite- 'time' limit of the correlator, in 
the sense that it is constrained by (95)-(96). The main 
assumption in the scaling form (93) is actually that the 
curvature of the correlator R{t, y) at the top of its cen- 
tral peak happens on a scale ^/y^ which corresponds 
to ^00 and is independent of 'time' and {c,D,T}. This 
assumption is not exact at all 'times', but we expect that 
it captures anyway the main features of the geometry of 
the correlator R{t, y) close to its central peak. 

Finally the three-point correlation function is assumed 
to scale in Dt and ^ in the same way as it naively does 
merely by counting the number of occurrences of rj in the 
definition (Dll) of R3{t,y) (i.e. inferred from R^rjfj, 
R3 ~ WnV aiid rescaling also the derivative dy): 
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with C4 = i?3(i,0) 



(97) 



Here C4 is also assumed to be a numerical constant. If 
this form is again not expected to be exact at all 'times', 
it can still be thought as a reasonable approximation 
provided the three-point correlator R'^{t,y) is analytic 
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FIG. 18. (Color online) Dt as a function of t {left: linear 
scale, right: logarithmic scale). In green (bottom curve), 
the solution of the full differential equation (99); in red (top 
curve), the solution of the equation without the non-linear 
KPZ term; the dashed horizontal lines are the corresponding 
large-'time' asymptotics; in purple dashed, the short-'time' 
asymptotics (102). Chosen parameters: c = D = T = £^ = l, 

Co = C2 = Cl = C3 = |, C4 = |. 

around y = 0. This last assumption is justified for in- 
stance in view of the zero-temperature and infinite- 'time' 
limit of the 'fiow' equation (34) (under the stationarity 
condition dtR — 0), which yields that the three-point cor- 
relator linit_>.oo R'^iiiV) is merely proportional to R'^[y), 
which is analytic around y = Q. 

Note finally that the scalings (92), (93), (97) are all 
compatible with the following rescaling in distribution 
(also inferred from rjfj ^ R) 

riit^y)= (98) 

for all values of the rescaling parameter a. A way to 
reformulate the definitions (94) and (97) of the numer- 
ical constants {02,03,04} is thus to identify those con- 
stants with the corresponding derivatives of the correla- 
tor R{t,y) taken at a ~ ^ and Dt = 1, which ensures 
their independence with respect to the other parameters. 



where = {£.cDY^^ is the same characteristic tem- 

perature as the one obtained in Sec. Ill C by scaling. 

Noting that C2 is defined in (94) precisely so that it ab- 
sorbs all the quantitative contribution of the geometry of 
R, the condition (96) guarantees that lini^_j.o / = 1- The 
^ — > solution of (100) is then / = so one obtains 

eventually that ~ 1 . This is valid for all values of the 

^ C2C3 

parameters if the {ci} are indeed parameter-independent. 
However those numerical constants are constrained only 
by the geometry of R{t, y), and we know from the numer- 
ical study in Sec. IV C that at large 'times' the correlator 
saturates to a -Rsat(y) ~ D^TZ{y) with only a slight T- 
dependence of the function TZ. 

Disregarding this possible but small modification of the 
numerical constant depending on T. the equation for the 
interpolating parameter finally reads 

Strikingly, it takes a form very similar to the equa- 
tions (46) and (70) obtained in Sec. Ill from the GVM 
approach, with an exponent 7 = 6 instead of 7 = |. In 
fact the value of this exponent only infiuences the specific 
monotonous crossover from the high-T regime, where 
the KPZ term has little infiuence (/ < 1) to the low-T 

asymptotics where / is linear in T as / ^ (^^) ^ T^- 
The value of 7 modifies the numerical constants in 
this last regime but does not infiuence the power-law de- 
pendence in the physical parameters, gathered in Tc(^). 
The equation (101) is a consistency check with respect to 
the two GVM predictions (46) and (70), and the strictly 
monotonous behavior of i?(i, 0) and D^o observed numer- 
ically in Sec. IV C. 



2. Evolution of Dt and prediction for D 00 and f{T,^) 

Substituting the rescahngs (92), (93), (97) into (91) 
transforms the equation for i?(t, 0) into an effective closed 
evolution equation for the amplitude: 



dtDt + -Dt = 



T ~ 



C4 1 
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(99) 



The non-linear KPZ term of the equation of evolution for 
F(t,y) corresponds to the term cx D^^'^ . The numerical 
solution of this equation is plotted in Fig. 18. 

To tackle the infinite- 'time' case, we have introduced 
in Sec. lllB-lllC the interpolating parameter / and 
identified it as the fuU-RSB cutoff in the GVM pre- 
dictions (cf. Appendix A). By cancelling the right-hand 
side of (99) and introducing the interpolating parame- 



ter /(c,D,T,C) 
equation for /: 



Doa/ ^ as in (45), one obtains a new 
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/) (100) 



3. Short-'time' evolution of Dt and saturation at tsat 

Leaving the infinite-'time' case, we consider now the 
opposite regime of short 'times' where the evolution (99) 
of Dt can be solved first in the^ absence of the non-linear 
KPZ term, predicting that Dt is linear in t at short 
'times': 



~ (t^o) coci D 

* — ^ — C2^ 
3c2 4 



(102) 



This behavior can also be checked by naked eye directly 
on (99) searching for a solution Dt oc t. Note the factor 3 
in the denominator, due to the two terms in the left-hand 
side of (99). Assuming that the solution Dt is also lin- 
ear at short 'times' while keeping the KPZ term actually 
yields the same result (since then D^^'^ <^ Dt <^ t°). 

One checks that this self-consistent hypothesis is cor- 
rect by solving (99) numerically (cf. Fig. 18). The in- 
dependence on T in the short-'time' regime (102) is also 
compatible with the direct simulation of the continuum 
DP (in a logarithmic scale of Fig. 9). The behavior (102) 



26 



thus holds in generality, and allows to define a saturation 
scale tsat at which Dt reaches its asymptotic value 



3c2 



tsat^^^f (103) 

CqCi T 



Accordingly to this equation, the saturation occurs ear- 
lier at higher temperatures, the thermal fluctuations ac- 
tually smoothing the evolution of the disorder correlator. 

Another consequence of the short-'time' behavior (102) 
deals with the short-'time' dynamics of the mean-value 
—2cFY{t,y) discussed in Sec. IV D and depicted in 
Fig. 13. Inserting the short-'time' Dt of (102) into the 
assumed scaling R{t,0) — of (93) and into the ex- 
act relation (33), one obtains from the initial condition 

Fvit,y) = that -2cdtFv{t,y) ^*=°^ hence the 

prediction 



VI. LINK TO EXPERIMENTS 

We discuss in this last section the consequences of the 
finite ^ > and its associated low-T regime for two spe- 
cific experiments. 

On one hand domain walls in ultrathin magnetic 
films^""^ are well described by the DES model of a ID 
interface defined in Sec. II A, and both their static geo- 
metrical properties and their quasistatic dynamical prop- 
erties (in the so-called 'creep' regime) are thus captured 
by the DP endpoint y{t) fiuctuations. 

On the other hand, a second instance of experiments 
encompassed by the KPZ theory is provided by interfaces 
in liquid crystals^^^^^, whose geometrical fiuctuations are 
directly described by the DP free-energy F{t, y) properly 
recentered. 



A. Temperature-dependence of the asymptotic 
roughness 



2c^V(^^*^°^^^t^ (104) 



This quadratic behaviour in t is compatible with the 
superlinear short-'time' regime observed in logarithmic 
scale in the inset of Fig. 13. However our numerical data 
do not clearly display such a behavior (delineated with 
dashed lines) on the whole range of temperatures, an ef- 
fect that we attribute to the thermal initial condition 
taken at Iq = 0.1 which alters the short-'time' power-law 
behavior that should be observed with the ideal initial 
condition iV(0, y) = at io = 0. 

This saturation 'time' is different from the character- 
istic 'time' t* — 



cD^ 



recently discussed in Ref. 65 for 
the evolution of the free-energy fiuctuations in the high- 
temperature regime and which corresponds to the Larkin 
length appearing by scaling at high-T or ^ = as dis- 
cussed in (54). This scale t* allows to draw apart a short- 
'time' diffusive and the large-'time' KPZ regime in the 
evolution of the fluctuations of F{t,y) (see also Ref. 46 
for a related study) , while the scale igat we examine in this 
section, singular in the limit — > 0, captures the short- 
'time' effects inherently due to the flniteness of ^ via the 
KPZ non-linearity. 

To summarize, at finite ^ the non- linear KPZ term 
does not modify the short-'time' regime but induces a 
saturation of Dt at shorter 'times' with increasing T 



and to an asymptotic value D^ < 



cD 



The high- and 
and 



low-temperature asymptotic regimes Doo — ^ 

— " ^S- are both independently well-controlled 

J- c 



D 



and the evolution equation (99) we presented thus allows 
to tackle the crossover from one regime to the other, but 
only in an effective way. An interesting open question is 
to provide a proper analytical derivation of the full tem- 
perature crossover and to fix the value of its exponent 7. 



As emphasized in Sec. IIIC and IV E, a prominent 
feature of the DP in a correlated disorder is that the 
amplitude of the roughness B{t) is modified by the mi- 
croscopic length ^ even at very large lengthscales (in the 
RM regime), provided that the temperature is lower than 
the characteristic temperature Tc = (^cD)^/'^. Interfaces 
in ferromagnetic thin films are a prototype system for 
the experimental study of ID interfaces, and in par- 
ticular their roughness exponent has been measured^'^ 
to be Crm ~ 0.66 in agreement with the KPZ exponent 
Crm = §• As estimated in Ref. 34 (Sec. VII B), the or- 
der of magnitude of Tc could be of room-temperature for 
these systems, which makes it even more relevant to de- 
termine whether they lie in the low-T or in the high-T 
regime. 

One could in principle distinguish between those two 
temperature regimes through the prefactor Arm of the 
asymptotic roughness at large lengthscales: 
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(105) 



as given by scahng arguments (66)-(67), predicted by 
GVM in (A15)-(A17), and consistent with the numeri- 
cal study of the roughness in Sec. IV E. 

The main problem regarding such a study is that the 
elastic constant c and possibly the disorder strength D 
may depend themselves on temperature, making it diffi- 
cult to characterize the low-T regime by a temperature- 
independent Arm, or to interpret a measure of the scal- 
ing in temperature Arm ~ T^^ with \i the thorn expo- 
nent. In any case, a change of regime in the temperature- 
dependence of Arm would provide a strong evidence for a 
low-T to high-T crossover. Promising experiments have 
actually been performed regarding the temperature de- 
pendence in ultrathin Pt/Co/Pt films, as analyzed in 
Ref. 69 where a fine study is devoted to the thermal 
rounding at the depinning transition. 
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In case of an elastic constant depending linearly on the 
temperature c = kT, we would expect a roughness am- 
plitude similar to obtained under the usual math- 
ematicians convention c = T, as exposed in Appendix G 
and with an additional dependence in k. 



B. Quasistatic creep regime 

The creep motion of ID interfaces, describing the qua- 
sistatic but non-linear response of the interface to an ex- 
ternal driving field, could also be an interesting bench- 
mark for the ^ > DP model predictions. 

Even though its description'^^''^'^ is not directly covered 
by the equilibrium statistical properties of the interface 
we have presented here, it happens that the characteris- 
tic lengthscales governing its scaling are actually believed 
to be the static ones. As detailed in Ref. 34 (Sec. VII B), 
those lengthscales are modified in the low-T regime, im- 
plying that the characteristic free-energy barriers scale 
differently with the temperature above or below Tc . This 
behavior could provide an experimental criterion to dis- 
tinguish between the low-T and high-T regime, especially 
since the exponent of the creep law has been successfully 
tested on domain walls in ultrathin magnetic films on 
several order of magnitude in the velocity 

A challenging situation would be that of an interface 
moving in a gradient of temperature, as studied numer- 
ically in Ref. 70, with a gradient spanning Tc itself and 
with an elastic constant c (or the disorder strength D) 
which could again depend or not on the temperature. 



C. High- velocity regime in liquid crystals 

Another experimental system where our approach 
might prove instructive is that of growing interfaces in 
liquid crystal turbulence^^"^^. The corresponding setup 
consists in a thin layer of liquid crystal subjected to a con- 
stant voltage U and to an alternating electric field which 
can generate two distinct turbulent modes (called 'dy- 
namical scattering modes') DSMl and DSM2, the later 
being more stable than the former. Starting from a dot 
(respectively a line) of DSM2 in a DSMl background, one 
thus observes the growth of circular (respectively flat) 
interface. We refer the reader to Ref. 68 for a complete 
account of the phenomena at hand. 

The fluctuations of this interface are actually remark- 
ably well described by the KPZ theory, providing a 
benchmark for its predictions not only about scaling ex- 
ponents but also about scaling functions. Contrarily to 
the case of magnetic interfaces, the fluctuations of the 
interface position are described by the random variable 
F{t, y) that plays the role of the disorder free-energy in 
the context of the DP, as defined by (12)-(13). Besides, 
for the liquid crystal interface, t is the physical time and 
y the longitudinal direction of the interface. 



For the circular interface, it has been shown that the 
experimentally measured C{t,y) defined by (19), once 
properly rescaled, matches very well the prediction of the 
KPZ theory at ^ = 0: C{t, y) can be fitted by the corre- 
sponding Airy2 correlator^^ . Although the lengthscale ^ 
of the disorder correlations is below the optical resolution 
of the experiment - no rounding of the cusp of C(t, y) is 
observed - it is still finite and should be relevant by in- 
ducing a crossover between two 'temperature' regimes to 
identify. The (non-centered) displacement h{t, y) of the 
interface evolves according to 

dth{t, y) - Woo + i^dlHt, y) + ^ [dyh{t, y)] ' + V{t, y) 

(106) 

where v^o represents the mean displacement velocity. 
Comparing this equation to the fiow of Fv{t, y) (25), one 
thus reads the correspondence A = ^, i/ = or equiva- 
lently c= \, T = It can be argued^'^'*^^ that in the 
experimental configuration of the growing dot, the pa- 
rameter A is directly given by the mean velocity A = Vao ■ 
This velocity itself is well described by a growing affine 
function of the applied constant voltage U , in the probed 
voltage range 26 V < C/ < 30 V (see Fig. 20 in Ref. 68). 
One may surmise that the parameter v, which describes 
the diffusive fluctuations of the interface, is independent 
of the constant voltage U (although this point is not dis- 
cussed in Ref. 68), and similarly for the amplitude D of 
the disorder fluctuations. 

If those assumptions are true, one can deduce from the 
correspondence c — T = ^ that the low-temperature 
regime of our description (where ^ matters in the scal- 
ing of the fluctuations amplitude) corresponds to a high- 
velocity A ~ Woo regime for the liquid crystal interface 
- and inversely for the low- velocity range being described 
by our high-temperature regime. A natural question is 
thus to determine whether the experiments are done in 
the high- or low-velocity regime. A test observable is 
provided by the amplitude of the fluctuations of h(t,y), 
denoted F in Ref. 68, and defined in the large-time regime 
from the scaling 

h{t,y)^^ v^t+{Tt)-^X2{y) (107) 
where X2 is the Airy2 process. In our notations, one has 

F = — (108) 
c 

This is read for instance from the factor which 
rescales the disorder free-energy in (73) with the choice 
a = {b / c^)^ ti oi (77). Equivalently, one has (Ft) 4 = E 
where E is the factor in (77) which rescales the free- 
energy in (78). From our intuition based on scaling ar- 
guments, the relation (108) should thus remain valid on 
the whole temperature range. The high-T and low-T re- 
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(109) 
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The crossover between those two regimes occurs at a 
characteristic Ac which is similar to the inverse of the 
characteristic temperature = ^/^D in the conventions 
of mathematicians (see Appendix G), as expected since 
T (X J and c — j. Strikingly, the crossover occurs be- 
tween a low-A regime (our high-T) where F is an increas- 
ing function of A and a high-A regime (our low-T) where T 
is a decreasing function of A. Note that the limit A — >■ 
in (106) is exactly the 'linearized' KPZ problem whose 
evolution has been solved at all 'times' in (40)-(41) and 
should thus be described exactly the microscopic disorder 
correlator R^{y) at asymptotically large 'times' (39). 

In Ref. 68, the authors have measured F as a func- 
tion of the voltage U and found a decreasing dependence, 
attributing it to a dependence in U of the parameters 
D or in the expression F'™ ^ = ^^-^ (tackled with 
^ = in our language, when C ^ can be neglected). In 
that spirit, (109) was first announced at the very con- 
clusion of Ref. 46 but with respect to the parameter v 
and with the two regimes being separated by the charac- 
teristic value Vc = {£,\^DY/^. We propose here another 
scenario where D and v are independent of U , and F fol- 

4 

lows its high-A expression F'^'shA _ d^^ ^ (confirmed 

A 3 ^ 3 

[e.g. by an independent measure of the parameter v), 
this would provide a clear evidence that the system lies 
deep in the high-A regime, and would to our knowledge 
constitute the first example of a phenomenon depicted 
by the 'low-temperature' KPZ regime. 



VII. CONCLUSION 

We have studied in this paper both analytically and nu- 
merically the consequences of a finite disorder correlation 
length ^ on the static properties of a ID interface depend- 
ing on the lengthscale - or equivalently a growing 1+1 
directed polymer at a fixed 'time' - in a random-bond 
quenched disorder accounting for a weak collective pin- 
ning. The two-point correlator R{t^ y) of the derivative of 
the disorder free-energy (the 'random phase' riv{t,y)) at 
fixed 'time' t emerged as a central quantity in the deter- 
mination of the dependence on the temperature T and 
the other parameters of the model (the disorder corre- 
lation length ^, the elastic constant c and the disorder 
strength D), as summarized in Sec. VA. 

The characteristic temperature Tc — (^cZ3)^/'^ sepa- 
rates two temperature-regimes in the 'time '-evolution of 
R{t,y), which is characterized by its amplitude Dt and 
the shape of its central peak TZ{y), from the point of view 



of our DP toymodel which focuses on the small trans- 
verse displacements \y\ (which are the most probable to 
be visited by the polymer and thus the more relevant 
ones). Although a priori the full shape of the correlator 
R{t,y) should matter, we showed that the asymptotic 
height of its central peak, parametrized by Z?oo(r, ^)/^ 
with ^ ~ ^, captures the main features of the crossover 
between low-T and high-T in the large-'time' random- 
manifold regime; the amplitude of the geometrical fluctu- 
ations in this asymptotic regime are actually not affected 
by ^ at high-T, whereas on the contrary at low-T it plays 
a crucial role and keeps this amplitude bounded. We 
showed that, surprisingly, the ratio /(T, ^) between the 
actual value of Dryo{T,S,) at finite ^ and its value ^ at 
^ = is directly related to a replica-symmetry-breaking 
cutoff parameter appearing in the GVM replica ap- 
proach, endowing this cutoff with an unsuspected phys- 
ical meaning (see Sec. IIIB). Moreover, Dt evolves at 
finite t according to two 'time' regimes, an initial regime 
and a saturation regime separated by a single 'time'- 
scale tsat which allows to understand how, depending on 
the temperature, the geometrical fluctuations depicted 
by the roughness B{t) display in turn two (at high T) 
or three (at low T) 'time' regimes. Finally the function 
7?.(?/) for 'times' above isat is strongly reminiscent of the 
microscopic disorder correlator R^(y) at high-T, but still 
needs to be characterized at low-T with some additional 
kernel to be determined analytically. 

We supported this picture through extensive numeri- 
cal simulations of the ^ > KPZ equation, which per- 
mitted in particular to evaluate Dao{T) numerically at 
a fixed value of f . In agreement with an effective equa- 
tion that we put forward for D^oiT), it appears that this 
quantity presents a crossover from low- to high-T, and 
not a phase transition that would be characterized by a 
non-analytic dependence of D^o in T. We cannot strictly 
exclude nonetheless the scenario of a phase transition be- 
cause of inherent numerical imperfections in the simula- 
tions and because the analytical equation on Drx,{T) is 
either effective (101) or obtained in the GVM approxi- 
mation (46) and (70). A rigorous procedure is needed 
to elucidate this alternative. Possible approaches encom- 
pass from an extension to the ^ > case of recent re- 
sultsi**'^'*"'^^^" in the ^ = KPZ class, to an adaptation 
of FRG arguments'^ ^"''•^ to the study of the correlator 
R{t,y). 

More broadly, it would be interesting to identify pos- 
sible connections between the 'time' equation of evolu- 
tion (34) of the correlator R{t, y), and the FRG flow equa- 
tion of evolution with respect to scale i of the renormal- 
ized disorder correlator in a FRG approach. Although 
the FRG framework holds perturbatively in e = 4 — d 
(hence e = 3 in our settings) the two varieties of flow 
equations bear striking resemblance. In particular, the 
high-temperature regime has been studied in Ref. 31 
by neglecting the non-linear contributions to the FRG 
flow, allowing to recover the same scalings as ours for 
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T ^ Tc- On the other hand, neglecting the non-linearity 
of the flow in our settings (Appendix E) indeed yields 
the correct high-temperature scaling Dao = ^ but does 
not provide for instance the correct roughness exponent 
Crm = I (when we consider the scaling of the negative 
bumps at large y, as in (41)). One would thus need a con- 
trolled expansion in order to draw a meaningful unified 
picture. Another related question deals with the T = 
singularity of the FRG renormalized correlator which ap- 
pears at a finite scale £c and plays an important role in 
the physical description of metastability in random man- 
ifolds. Our result should help to understand in which 
order the limits T — and — are to be taken. Last, 
a FRG approach, through the possible existence of a zero 
and a finite temperature fixed point, may allow to dis- 
tinguish between the crossover and the phase-transition 
scenari. 

The richness of the KPZ universality class also allows 
to translate the occurrence of low- vs high-temperature 
regimes into different languages. We have illustrated this 
fact in the analysis of a liquid crystal experiment^^"^^ 
(Sec. VI) where the low-temperature regime corresponds 
to a high-velocity regime - yet to be ascertained exper- 
imentally. On the mathematical side, keeping ^ finite 
amounts to generalize the Airy2 process in a non-trivial 
way (see also Ref. 46). In the language of replicae^^, 
it corresponds to solving the problems of bosons with 
attractive but non-S interactions in one dimension, for 
which the ^ = Bethe Ansatz solution is not known to 
generalize. The same search for an extension also applies 
to the Airyi process, which describes a point-to-line DP 
problem (with 'fiat' initial conditions), which finds an 
experimental incarnation for instance in fiat interfaces in 
liquid crystals. 

One may finally wonder how the existence of a low- 
temperature regime extends to phenomena which do not 
fall a priori into the KPZ universality class. A basic 
assumption made at the very start is that the disorder 
is Dirac 5-correlated along the longitudinal 'time' direc- 
tion. The disorder is nevertheless always correlated in 
both directions in physical systems; one generalization of 
our results would be to take this property into account. 
Such correlations may also account for overhangs present 
in interfaces, once smoothed out by a change of scale. 
The generalization is non-trivial in the sense that the 
equations of evolution become non-local in 'time'. More 
broadly, the cases of non-RB disorder and/or higher di- 
mensions, where a mapping to a directed-path problem 
is not always possible, are still open. 
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Appendix A: Reminder of previous GVM roughness 
predictions 

For completeness, here we recall and adapt the rough- 
ness predictions for the static ID interface obtained in 
Ref. 34, first assuming that the disorder correlator R^{x) 
in (1) was a normalized Gaussian function of variance 
2^^ , then using the replica trick in order to average over 
disorder and finally performing a Gaussian Variational 
Method (GVM) with a full replica-symmetry-breaking 
(RSB) variational Ansatz as introduced by Mezard and 
Parisi in Ref. 45 and 71, further investigated also by 
Goldschmidt and Blum in Ref. 43. We emphasize in 
particular the role of the interpolating parameter f{T, ^) 
between the high- and low-temperature regimes. 

For the full DES model of the ID interface, we had 
obtained for the variance of the relative displacements 
B{r) = (Au(r)2) at the lengthscale r: 
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Sdis(r) 
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where the four DES parameters {c, DjT, ^} are respec- 
tively the elastic constant c (elastic energy per unit of 
length along the interface), the disorder strength D (the 
typical amplitude of the random potential) , the tempera- 
ture T and the disorder correlation length ^ (or width of 
the interface), tq is the Larkin length introduced in (58) 
and marking the beginning of the asymptotic 'random- 
manifold' regime, Vc{T,S_) the fuU-RSB cutoff and Tc the 
characteristic temperature separating the low- and high- 
temperature regimes. 

As for our DP 'rounded' toymodel, assuming that the 
effective disorder correlator of ryy (t, y) (20) is of the form 
R{t, y) — D ■ TZ^iy) we have performed a GVM procedure 
on the following statistical average obtained from (22) 
using replicae to average over disorder: 
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With the function TZAy) the same normalized Gaussian 
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function of variance 2^^, we have obtained for the vari- 
ance of the DP's endpoint fluctuations -Bdp(^) = (2/(i)^) 
after a growing 'time' t: 



BDp{t < tc) = — + 



tr. = 
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where, as before, the four parameters |c, T, are 
respectively the elastic constant c, the temperature T, 
the effective strength of disorder D, the effective interface 
width and Mc(T,0 the fuU-RSB cutoff. 

However, the second GVM computation was performed 
at a fixed 'time ' t, so the effective parameters Dt and 
have a priori also a 'time' dependence. Assuming that 
they both saturate quite quickly, we can safely compare 
the two sets of predictions (A1)-(A5) and (A7)-(A10) 
with the translation of 'time' t into the lengthscale r, the 
approximation ^ ~ ^, and the identiflcation of tc with tq 
for the Larkin length. As for the full-RSB cutoffs, the def- 
initions f'^iT,^) = |z;,(T,0 and f°HT,0 = |wc(T,0 
yield then two similar equations for the interpolating pa- 
rameter /(T, ^), respectively (70) and (46), if we im- 
pose by hand D = ^/'"^(T,^) in order to match the 
Larkin lengths tc and tq. Their numerical discrepancy 
can be safely attributed to the GVM approximation. The 
structure cx {T/Tc)^{l — /) stems in both cases from 

the comparison of^thiTf = {^f and = it\,{Tc?, the 
equation (A4) being initially of the form: 
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As discussed in Ref. 1 1 , the Larkin length is a physical 
benchmark for the roughness, firstly as the beginning of 
its asymptotic 'random-manifold' regime, secondly with 
its relation to the maximum value of the GVM self-energy 
[a] {vc) and thus to the full-RSB cutoff itself Vc{T, ^), and 
thirdly consistent with its original definition by Larkin^^ 
as the lengthscale at which the typical relative displace- 
ment of the ID interface corresponds to its effective width 
B{Lc{T, $)) « fcff (r, if which fixes the ampHtude of the 
RM roughness B{r > Lc) « A{c, D, T, i)r^^^: 
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Using (A4), (A3) and (66), we give the GVM predictions 



for the quantities 
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Appendix B: Statistical Tilt Symmetry (STS) 

We justify in this appendix the decomposition of the 
free-energy Fy{t,y) defined by (12) into the sum of a 
disorder-independent term Fv=o{t,y) and a translation- 
ally invariant term Fv{t,y) (13-16), from the comple- 
mentary viewpoints of path-integrals and of stochastic 
differential equations. This symmetry arises from three 
ingredients: the precise form of the elastic energy den- 
sity ^{dtyf , the invariance in distribution of the dis- 
order V{t^ y) by translation along y and the continuum 
nature of the transverse direction y. We refer the reader 
to Ref.^®'^^"^'* for previous discussions of the STS. 

The weight of trajectories {y{t)} starting in (0,0) and 
arriving in (ti, yi) can be compared to the weight of those 
arriving in (ii, 0) from the change of coordinates: 



y{t)^y{t)-ft 



(Bl) 



with now the initial and final conditions y(0) — y{ti) ~ 0. 
Introducing the tilted disorder 



7^lV{t,y) = V{t,y+ft) 



(B2) 



and performing the change (Bl) in the path-integral of 
the unnormalized weight (3) one obtains: 

Wv{ti,yi) = e-^i/(2^-(*i»iy.^nv^(0,ti) (B3) 

thanks to the form ^{dtyf of the elastic energy den- 
sity which allows to identify a translated path integral 
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over {y{t)} and to single out the thermal contribution 
with Bt\i{t) = The measure of the path-integral re- 
mains unchanged {'Dy{t) = T?y{t)) since the change of 
variable (Bl) is a translation. Using the definitions (13), 
this equality writes for the disorder free-energy 



Fv{ti,yi) 



(B4) 



This means that the whole dependence in the arrival 
point 2/1 of the disorder free-energy can be absorbed 
into a tilt Ty^ of the disorder potential. Using that at 
fixed final time ti the disorder is translation-invariant 
in y i.e. P [V] ^ P [TylV] one obtains the symmetry 
P[Fv{t,y + Y)] ^P[Fv{t,y)] as announced in (17). 
This translation invariance in the y-direction is valid for 
any functional of Fv(t, x), and yields in particular for the 
/c-point correlators: 



Fv{h,yi+Y)...Fv{h,yk + Y) 



Fv{ti,yi)...Fv{ti,yk) 



(B5) 



The same result (17) can also be deduced from the 
strict point of view of stochastic differential equations. 
Note first that the initial condition .^^^"'^'i = 5{v) trans- 

lates for Fy as (0, y) = which is trivially translation- 
invariant along the y-direction. Showing the STS (17) 
thus amounts to checking that this property is preserved 
in time through the evolution equation evolution (26) for 
Fy. Consider for this purpose a 'Galilean transformation' 
of Fv{t,y), which consists in defining Fy(t,y) through 



Fvit,y) = F^{t,y-vt) 



(B6) 



where v represents the 'velocity' of the tilt y i-^ y + vt. 
One sees directly from (26) that the terms in v in the 
equation of evolution for Fy(t,y) compensate between 
left and right hand side: Fy{t, y) verifies the same equa- 
tion of evolution as Fv{t,y), but in a tilted disorder 
V'"{t,y) = V{t,y + vt). Since the initial condition is 
also the same, one obtains: 



FUti,yi)^Fv^{h,yi) 



(B7) 



Choosing v = yi/ti in (B6) and (B7) yields again (B4). 

Another but less general incarnation of the STS arises 
when considering the geometrical fiuctuations instead of 
the free-energy. Defining the 'generating function' of the 
moments of y{ti) 
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where the final-time condition is free (or equivalently, 
y{ti) is integrated over), the disorder average of the vari- 
ance of the endpoint y(ti) is given by: 



{yihr),^{y{hr)-{y{h)y^di 
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We note that performing an appropriate tilt encoding 
the thermal roughness, namely y(t) = y{t) + XBthit) one 
obtains for the generic function: 

W^ih) = e^"'^*^)^'/^ w^=°ih) (BIO) 

where V{t, y) = V{t, y+ ^t) is a tilted disorder. Taking 
the logarithm and averaging over disorder, one obtains 
for (B9): 



(y(iinc = ^-5th(ti) 



(Bll) 



and {y{h)k)^, ^d'^\^^^\ogW^{ti)^0 for fc > 2. It 
would be tempting to conclude from these averaged 
cumulants that the average of the whole distribu- 
tion Pv{ti,yi) is a normal law J\f{0, Bth{ti)); this is 
only true in the trivial case without disorder where 
Pv^o{t,y) = Pth(<,2/) = AA(O^Bth(ti)). Indeed the nor- 
malization at fixed disorder Wyiti) of definition (5) pre- 
vents a direct disorder average on (B3). As for the rough- 
ness B(t), the property (Bll) can be reformulated in an 
intrinsic way comparing the cumulants with respect to 
the thermal, disorder and the joint disorder-thermal dis- 
tributions: 



Bit)^{y{tr)^{ym^+{yit)y 



Bdis(t) = B(t)-Sth(t) = (2/(i))' 
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Physically this decomposition allows to focus on the pure 
disorder contribution of the roughness _Bdis(i), which is 
equivalent to the full roughness at large 'times'. Such a 
simple relation is possible for the second cumulant thanks 
to the following generic identity valid in presence of two 
probability laws pi and p2 (respectively thermal and dis- 
order distribution in our case): 



Varpjopi = o Var 
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Appendix C: Derivation of the Feynman-Kac 
'time'-evolution equations 

We rederive in this appendix the 'time'-evolution equa- 
tion (23) of the weight Zv(t,y) = defined with 

^ ' V V , t,/ Wv = o{t) 

respect to the pseudo free-energy Fv{t,y) in (12). Given 
in Ref. 29, this evolution equation is the starting point 
of the Feynman-Kac equations of Sec. II D which de- 
fine univocally the pseudo free-energy quantities Fv{t, y), 
Fv{t, y) and its random phase riv{t, y) with their ad hoc 
initial conditions. However, the normalization of the 
weight Zv(t,y) is not conserved in presence of disorder, 
so it requires a careful treatment in order to yield its 
continuous formulation (23). 

Using infinitesimal 'propagators' in a path-integral for- 
mulation of the weight Zy'^^ , we derive thereafter the 
evolution equation first using its propagation equation 
in continuous time and secondly constructing explicitly 
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the weight with the propagators in discretized time. We 
rederive this result for completeness, to pinpoint the re- 
quired hypotheses and the possible issues in a generahza- 
tion regarding the form of the elasticity. 



1. Propagation equation in continuous time 



We first define the following propagator (with /3 = 7^ ) : 



Zviti,yi\to,ya) = 
n[y,V;to,h] = 



y{ti)=yi 
y{to)=yo 



Vy{t) e-'3«fe'^;*o.til (ci) 



dt 



-{dtyf + V{t,y{t)) (C2) 



which represents the weight of trajectories starting in 
at time and ending in yi at time ti > tQ. Note that 
by definition Zy{t,y) = yjO, 0). Paths y{t) are 

weighted by a measure 'Dy{t) ensuring that for purely 
thermal paths the weight is normalized. Explicitly, with 
the normalization Wv{t) defined by (4), the measure 

-Dyit) 



-Dyit) 
in (CI) ensures that 



dyi Zv=oiti,yi\to,yQ) = 1 



(C3) 



(C4) 



a property which is not true anymore for any V. Here, 
/ dy denotes dy as in the rest of this appendix. The 
advantage of using this choice of normalization is that 
Zv{t, y) obeys the so-called stochastic heat equation (23), 
which can be shown using the Feynman-Kac formula. 

From its definition, we see that the propagator (CI) 
presents several useful properties. For ti — to it goes to 
a Dirac delta: 

lini Zv{ti,yi\tQ,yo) ^ 5{yi- ya) (C5) 

since the trajectory endpoint yi remains very close to its 
departure point at small times; to be more precise, we 
read from (CI) the expression of the infinitesimal prop- 
agator 5Zy, valid for ti close to tg- 
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The prefactor ensures the normalization condition (C4) 
and also yields (C5) in the limit ti — > to. It is rather 
important not to overlook this prefactor since it ensures 
that the infinitesimal propagator (C6) evolves in time ac- 
cording to the forward and backward equations (already 
close to the final one on Zy): 
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dtJZy = - 
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Last, the non-infinitesimal time evolution is described by 
the propagation equation 

Zv[ti,Vi\h,yo) = J d,yZv{ti,yi\t,y)Zy{t,y\to,yo) 

(C9) 

which expresses that the path integral (CI) can be cut at 
any time to < t < ti provided the intermediate values y 
of the path at time t are integrated upon. Using (C5), the 
limit t — > ti of the equation of propagation (C9) yields a 
trivial identity. To go further, we can differentiate (C9) 
with respect to time t and then use that for t close to ti 
the infinitesimal propagator verifies the backwards evo- 
lution (C8) in order to write 
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dy 
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dy SZy{ti,yi\t,y)dtZy{t,y\to,yo)} 
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Integrating by parts and taking the limit t ^ ti thanks 
to (C5) finally yields the expected stochastic heat equa- 
tion 

1 .2 



dtZv{t,y\to,yo) 
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Note that a propagator obeying such an equation of 
evolution cannot keep its normalization constant since in 
general 



dt J dyZy{t,y\to,yo) ^ ~^ J dyV{t,y)Zy{t,y\to,yo) 

(C12) 

is non-zero. We also remark that considering an elas- 
tic energy including higher powers of dty than {dtyY 
in (C2) would be problematic for finding the equation of 
evolution: the infinitesimal propagator (C6) would con- 
tain terms of the form making it not obvious 

to determine the equivalent of (C7-C8) and deriving the 
evolution corresponding to (Cll). 



2. Explicit propagator in discretized time 

A second approach to explicit the normalization prop- 
erties of the path integral is to work in discretized time. 
A path is going from j/i to yf between time t{ and if in 
N time steps 6t = so at times tk = + k^^^. We 

define the weight (or probability density) of a free path 
as 
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= n 9{St,yk+i-yk) (C13) 

0<fc<Ar 

where we have denoted the microscopic propagator by 



9{t,y) 



-/3cl 



(C14) 



33 



One defines the discrete equivalent in Ito's discretization 
to the continuous propagator (CI) as 

2v [tf, yi\U, Vi) ^ Jdyo... dyN Viyo ■ ■ ■ Vn] 

'x S{yQ-yi)S{yN -yi) (C15) 



X exp 



yk] 



0<k<N 



with the expectation that this result does not depend 
on A'' in the large N limit. The following decomposition 
makes the link with the continuum formulation (C1-C3): 
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At y = the normalization corresponding to (C4) still 
holds: the writing (C13) is a product of normalized prob- 
ability densities. The propagation equation (C9) is read- 
ily verified: for all intermediate times t = tg one has 

j dyZv{tf,y{;t,y)Zy{t,y;ti,yi) 

= J'^y Jdyo- - ■ dyedy'i . . . dyN S{yo - yi)S{ye - y)S{y'e - y) 
X S{yN -yf)P{yo---yN]exp 5tV{tk,yk) 

0<k<N 

= J dyo... dye... dyN S{yo - yi)S{yN ~ yt) 

X P[yo ■ • ■ 2/iv] exp - /? ^ StV{tk,yk) 
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The derivation of the stochastic heat equation can 
then be made explicit: fixing now y\ and (and skip- 
ping them in the all following weights Zy (. . . j/i)), 
one compares two histories ito,yo] . . .',tiy,yN) and 

(^0, 2/o; • • ■ ; tN: yN;tN+i,yN+i) 

2y^^itN+i,yt) = I dyN+i . . ■ dyoS{yo - yi)S{yN+i - yt) 



~l3St E V{tk,yk) 
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where we have isolated the contribution of the last 
time step. Expanding the exponential g-f^^tvitN+i,VN+i) 



(which is valid at minimal order in St), yields at order St: 
Z^+\tN+i,yi) - (tM, yi) ^ -Stpvz^ {tN,yi) 
dyN+i . . ■ dyoSiyo - yi)S{yN+i ~ yt) 



X [giSt, yi - yN) - S{yi - yjv)] 

xV[yo...yN]e °^'=<" (C19) 
where the third line reads 

9{yi-yN,St)-5{yf-yN) 

= g{St,yi-yN)-9{0,yf-yN) (C20) 
^5tdtg{5t,yi-yM) (C21) 

Using now dtg{5t, yt - yn) = ^c^ltdi^t, yf - yN) and, in 
the integral, dy^ — dy^^^-^ and integrating by parts we get 
finally 

Z^+\tN+i,yi)-Z^itN,yi) _ 
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which corresponds to the continuum equation (Cll) in 
the limit St — > 0. Note moreover that from (C18) the 
following exact recurrence equation for the propagator 
can be read: 



'{U + St,yf) = 



dy g{St, yi - y)e-''^*^(*«>^)Z^ (tf, y) (C23) 



and is the continuous analogue of the transfer matrix 
equation for a directed polymer constrained on a dis- 
cretized lattice. 

Note that similarly to the continuous case, the specific 
form of the short-range elasticity implies the Gaussian 
form of the microscopic propagator (C14), which thus 
satisfies the diffusion equation leading to the discrete 
stochastic heat equation (C22) . However a different elas- 
ticity will in general not be Gaussian and consequently 
radically change its evolution equation. 



Appendix D: 'Time'-evolution equations of averages 
using the functional Ito formula 

Aiming at the derivation the evolution equations for 
dtC{t,y) and dtR{t,y) given in Sec. HE, we first present 
in this appendix the functional Ito formula''^ applied to 
a field X{t, y) obeying a generic Langevin equation, then 
we particularize it to the case of multi-point correlators, 
obtaining finally the flow equations (34)-(35). 

The Feynman-Kac evolution equations (25), (26) and 
(27) of Fv, Fv and rjv take the form of a generic 
Langevin equation for a field X(t,y): 



dtXit,y) ^g[X{t,y);t,y]+Vit,y) 



(Dl) 
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where V(t, y) is a centered Gaussian noise with correla- 
tions 



V{t, y)V{t', y') = D5{t' - t)n{y' - y) (D2) 

For X{t, y) being y), Fv{t, y) and r?y (t, y) one reads 
g respectively as g[F] from (F2), G[F] t, y] from (F4), and 
t, y] from (F6). For F and F one has 7^(?/) = R^{y) 
while 7?.(j/) = —R'^{y) for X — -q (note that this par- 
ticular notation is specific only to this appendix). In 
what follows one denotes y] for short instead 

of Q[X{t,y)]t,y\. Our aim is to deduce, from the 
Langevin equation (Dl), evolution equations for the sta- 
tistical average of functions of {X{t,yi)} at different 
points {j/i}; such as the two-point correlation function 

X{t,y^)X{t,y2). 

This is fairly straightforward to obtain the time- 
derivative of X(t,y) by directly averaging (Dl) which 
yields dtX{t,y) = Q[X\t,y] {e.g. (32) and (33)). The 
same cannot be applied to X{t,y)'^ since, when writing 



dtX{t,y)^ ^2Xit,y)dtXiLy) 



(D3) 



= 2X{t, y)g[X; t, y] + 2X{t, y)V{t, y) (D4) 

one cannot easily eliminate V from the last term. To 
tackle such correlation functions, one may use the 
functional Ito formula which reads as follows, for the 
Langevin equation (Dl) with continuous argument y and 
non Dirac delta correlated random potential V: 



dtg[X] = Jdyg[X;t,y] 



6X{y) 



+ f/«W-.)^^^g|L (D5) 

where g[X] is a functional of X. For g[X] = 0{X{yi)) 
where yi is fixed {e.g. an observable depending on the 
sole DP endpoint), one thus finds 



dtO{X{t, yi)) - g[X; t, y,]dxO{X{t, y,)) 



D, 



+ -n{0)dj,O{X{t,y,)) (D6) 



For 0{X) ^ X^ one finds 



dtX{L yi)2 = 2g[X; t, yi]X{t, yi) + 071(0) (D7) 

which is the correct form of (D4). Note that the result 
is singular for ^ — ^ in our cases of interest TZ = and 
n = -R'^. 

Another example is provided by the computation of 
the time evolution of the average of multiple-point cor- 
relators, for which (D5) yields: 



dtX{t,yi)Xit,y2) = g[X;t,y,]X{t,y2) 



(D8) 



g[X;t,y2]X{t,yi) + Dn{y2-yi) 



which yields back (D7) for yi = y2- More generally, one 
has, noting diO (respectively 820) the derivative of O 
with respect to its first (respectively second) argument: 



dtO{X{t,y,),Xit,y2)) 



-g[X;t,yi]diOi. 
D 
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.)+g[X;t,y2]d20{...) 
27^(y2-yl)5l2 + 7^(0)922 



0(...) 

(D9) 



We now derive the evolution equation for the corre- 
lator R{t,y) of riv{t,y). We first explicit some useful 
parity symmetry. The equation for riv{t,y) is the same 
as for —r]yR{t,—y) with a reflected disorder V^{t,y) — 
V{t, — y). This proves that 

ilv{t,y) = -7jvn{t,-y) (DIO) 

at all times. Since the distributions of V and are 
the same, one can replace in averages every 'i]{t, —y) by 
—ri{t^ y) without changing the result (so from on we skip 
the index V). In other words, 7]{t, y) is an odd function of 
y in distribution. We now define a three-point correlation 
function 



R:i{t,y)^ri{t,y)Mt,^) 



(Dll) 



To simplify the notations, and since one only considers 
one-time observables at time t, we now drop the depen- 
dence in t and denote the derivation with respect to y 
by a prime. Using the noted parity (which also extend 
to the derivatives; e.g.: ri'(y) is an even function of y 
in distribution^^) together with the statistical invariance 
by translation, one finds for instance, by translating all 
arguments by — y 



R'\y) = v''iyHof^'v''io)rihyr^ 'ri''{oUy) (D12) 



(pa 



2^3(2/) 



(tr.) 



V'ivHyHO) = r?'(0)77(0)r;(-y) 



(pa 



77'(0)77(0)?7(y) 



(D13) 



We are now ready to determine the time-evolution of 
R(t,y) combining (D8) and (F6): 



dtR{t,y)^dtv{0)Tj{y) 
(m)T 
2c 



(D14) 



v"i0Hy) + v"iyH0) 



v'ioMOHy) + v'iyMy)v{0) 
1 



t 



277(0)r/(y)+y77'(y)?7(0) ~ DR'Hy) 



We eventually recognize thanks to (D12) and (D13) that 
T 



dtR{t,y)=-dt.R{t,y) - -dyR3{t,y) 



R{t,y) + dy{yR{t,y)) -DR'Uy) 



(D15) 
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This equation is valid at all times and would in principle 
allow to studying of the 'flow' of R{t, y) starting from its 
initial condition ^(0, y) = 0. Due to the non-linear KPZ 
term, it is however non-closed on the two-point correla- 
tion function y) and brings into the game a three- 
point correlation function Rz{t, y), making it necessary to 
solve the full hierarchy of equations for the n-point func- 
tions to determine R(t, y). Yet, using a scaling Ansatz in 
?/ = 0, the equation (D15) still enables to determine the 
time evolution of the height of the R{t^ y) in y = (see 
Sec. V). 

Similarly, defining the three-point correlation function 
for F 



C3{t,y) = -2[F{t,y)^F{t,0)] [F'{t,0)Y (D16) 

one obtains the flow of C{t,y) using (F4) in (D9) 
with _0{Xi,X2) = (Xi - Xif, = F{t,y) and 

X2 = ^(i,0): 



d,C[t, y)^- [C"(t, y) ~ C"{t, 0)] - y) 
--Ci[t,y)-2D [R^{y)-R^m 



(D17) 



Appendix E: Solution of the linearized dynamics of 
F for a generic disorder correlator R^iy) 

In this appendix, we determine the explicit form of 
the correlator C{t, y) (resp. R{t, y)) of the disorder free- 
energy F (resp. of the random phase 77), in the ap- 
proximation where the 'time'-evolution equation of those 
quantities is linearized. The crossover from flnite to infl- 
nite 'time' regime is discussed. 

The evolution equations (35) for C{t, y) and (34) for 
R{t, y) are not closed because of the three-point corre- 
lation functions C^{t,y) and R3{t,y). It is yet instruc- 
tive to solve those equations in the approximation where 
those three-point functions are set to zero. An equiva- 
lent alternative approach is to solve directly the equa- 
tion (26) for F or (27) for 77 by neglecting again the 
non-linear terms in those equations (see also Appendix C 
'Short-time dynamics (diffusive scaling)' of Ref. 46, for 
the explicit case of Gaussian function for the microscopic 
disorder correlator R^{y)). We denote by C'''"(<, y) and 
^''"(t, y) = ^^^(^''"(t, y) their solutions, which are ex- 
pected to be valid either at small 'times' for all y (be- 
cause the initial condition ensures those functions vanish 
uniformly at 'time' 0) or at all 'times' but small y in the 
high-temperature regime (see Sec. VB for a discussion). 

The equations at hand are linear and are thus solved 
e.g. using Green functions, and the solution takes the 
form 



(El) 



C'"^(t,y) / dw Kt{y,w)R^{w) 




/l.O 



FIG. 19. Kernel k{y, w) defined in (E5) linking the Hn- 
earized free-energy correlator C''"(t,j/) and the disorder cor- 
relator with the combination of (E3)-(E4). 



where the kernel reads Kt{y, w) — dsKs{y, w; t) with 
Pet 1 



Ksiy,w;t) = 



2^7^ y/l3cst{t - s) 



(E2) 



■ get /jc itw-syy 



I3c (tm + ay)^ 
g 4 ts(t-s) 



The scaling analysis of this expression may be obtained 
by setting s — tr. To this end, we assume the natural 
rescaling R^{ay) = a~^R^/a{y) of the disorder correla- 
tor. One obtains that the free-energy correlator rescales 
purely diffusively as 

C'^{t,y) = ^ C^(^^\ (E3) 

with as usual B^_]^{t) — ^ and with the scaling function 

C-^{y)^ j dwk{y,w)R-^(w) (E4) 
dr 1 



K{y,w) 



y/T{l - t) 



(E5) 



2e 4t(1-t) _ g 4t(1-t) _ g 4x(1-t) 



Under this scaling, the whole {c, D, T}-dependence is ab- 
sorbed in the prefactor The scaling kernel K{y,w), 
illustrated on Fig. 19 is continuous but non-analytical on 
the lines |y| = |u;|, as can be seen from the direct compu- 
tation of (E5), which, using the symmetry by even parity 
is expressed for y > as: 



(y>0) 



dw K"" (y, w)Ri;{w) 

roc 

+ / dw jk^{y,w)R^(w) 
■ly 



^-^^ = 1 - Erf - ^ ^2 - Erf ^ - Erf 
2^ ^ - - 



(E6) 



f' 



l_Erf.^-^(2 + Erf^^-Erf 



20F 



36 




-10 -5 5 10 -5 5 



FIG. 20. (Color online) The finite-'time' correlator R^'"^it,y) 
(thin purple lines) for the 'Gauss' disorder correlator R^{y) = 
Rf^^^^{y), plotted as a function of y for different 'times', com- 
pared to its infinite-time limit R^{y) (thick red line), with 
^ = 1, c = 1, -D = 1, and T = 1. Larger 'times' correspond to 
lighter colors. Left inset: same behavior for (5''" (t, y). Right 
inset: same behavior for ^dy&'^{t,y) — J^^ dy' R^^^{t,y'). 

Those expressions describe through (E3) the complete 
transition from the initial regime where the correlator 
is close to zero to the infinite-time asymptotic regime, 
where one should recover limt^ao R^^'^{t,y) — ^R^{y) 
(see Sec. F3). 

As we now detail, this infinite- 'time' limit is however 
not obvious to extract from (E6) and is in fact directly 
related to the non-analyticity of K{y,w). Carefully inte- 
grating by part and differentiating with respect to y leads 
to the following expression, valid provided that R(^{y) is 
bounded at infinity 

Ic'iy) ^^i°'4-i)(y) + r dw we-''^^+y^R^f'\w) 
2 Jo 

/•OO 

- / dw we-'^^'^-y^ R^f^\w) (E7) 

Here i?^ ^\y) is the primitive of R^(jj) which vanishes 
in 0. It verifies the scaling relation _R^~^''(aj/) = i?^^j^'(?;) 
and its small ^ limit is half of the Heaviside step func- 
tion lim^_j.Q i?^~"^^(j/) = |6(y), provided now that R^{y) 
describes a RB disorder. Its occurrence as the first term 
of (E7) arises from the jump of the slope of K{y, w) in 
y = w, depicted in Fig. 19. Using those properties one 
obtains 

lim l^C ' ; ( —^=J==) > ^ dwwe-'^^0(iD) 

= R'f^\y) + ^lim ^J^dw we-'^^'^+y"^ R^~^\w) 
- j dw we''"^''-y'^ R^';^\w) I 

V ' 

^R^'\y) (E8) 



FIG. 21. Same correlators as in Fig. 20, now corresponding to 
the 'CubicS' disorder correlator Ri^{y) = Rf""""'"^ (y) , defined 
in (I2)-(I3). Parameters are ^ = 2, c = 1, D = 1, and T = 1. 



where we have denoted for short y* = y / ^Jli^^Jf) and 
^* = Bt\i{t). Differentiating with respect to y, one 
obtains the expected result limt_j.oo R^™{t, y) — ^R^{y), 
while keeping t finite yields the decomposition (40-41) of 
R{t,y) announced in Sec. Ill A. The short to large-'time' 
behavior of the correlators is illustrated in Fig. 20 for the 
Gaussian disorder correlator R^{y) = TZ^'^^^^iy) defined 
in (88). 

As depicted in Fig. 21, one may also measure numer- 
ically the correlator R}™{t,y) for the disorder correlator 
R^{y) = i?CubicS(y) defined in (I2)-(I3). As in Sec. IV B 
we fitted the average of this correlator over large 'times' 
at £^ — 2, c = 1, D ~ 1 and T = 1, using the 'Gauss' and 
the 'SincG' fitting functions. Parameters are found to be 

^slr^/^s^ar^^ ^ 0.89 ar^/Sr^^ ^ 3.27 (e9) 

which correspond to the same results as in the high- 
temperature regime of Sec. IV B (see Table I). We 
checked that this numerical results depend only slightly 
on the value of the temperature, decreasing with T. 



Appendix F: Fokker-PIanck equations for the pseudo 
free-energy 

Starting from the path integral formulation of the 
pseudo free-energy, the average over thermal fluctua- 
tions yields the so-called 'Feynman-Kac' equations for 
dtWv{t,y), dtFjj{t,y), dtriv{t,y)- In this section, we re- 
examine the possible steady-state solutions of the Fokker- 
PIanck equations, obtained after the disorder average 
over the random potential, and try to generalize them 
from the uncorrelated disorder (^ = 0) to the case of a 
random-bond correlated disorder (^ > and short-range 
correlator) . 
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1. FP equations for Fy, Fv and r]v 

The pseudo free-energy Fv{t, y) follows the KPZ equa- 
tion (25): 



dtFv{t,y)^g[Fv{t,y)] + V{t,y) 

g[F]^'^dlF{y)-^[dyF{y)f 



2c 



2c 



(Fl) 
(F2) 



and similarly the disorder free-energy Fv{t, y) follows the 
tilted KPZ equation (26) depending explicitly on {t,y): 



dtFvit.y) = g [Fv{t,yy,t,y] + V{t,y) 

y 



g\F;t,y =g\F 



dyFiy) 



(F3) 
(F4) 



as its corresponding random phase riy{t, y) = dyFyit, y): 
dtwit, y) = Gri [vv{t, y);t, y] + dyV{t, y) (F5) 



T 



dy[v{yf 

2c 



a, 



t 



riiy) (F6) 



Combining the Feynman-Kac equations (25)- (27) and 
the random-potential disorder correlator: 



V{t, y)V{t', y')=D- 5{t - t') ■ R^{y ~ y') 



(F7) 



the 'time '-evolution of the free-energy distribution 
V [F, t] is then given by the functional FP equation (ob- 
tained e.g. from Ito's lemma): 



dtV[F,t] 
D 



dyT^A~Q\F]T'\FA] 



5F{y) 
dydy' ■ R^{y - y' 



S^V[F,t] 
SF{y)6F{yO 



with those two terms stemming respectively from the de- 
terministic operator g and the remaining stochastic term 
in the Feynman-Kac equation. The distribution V [Py, i] 
follows the same functional FP equation, with the tilt of 
the deterministic operator g [F] i-^Q [F; t, y] , whereas 
the random-phase counterpart satisfies: 



dtV[v,t] - 

J dydy' -R'^iy-y' 



dy^^{-gn [mtMPliiA] 



D 
2" 



S'q{y)Sr]{y') 



(F9) 



The detailed notations Fy{t^y), Fy{t,y) and i]y(t,y) 
have been simplified to the random functions F{y), F{y) 
ad ri{y) of respective distributions V [F,t], V [Fjt] and 
V [t], t] at a fixed 'time' t, though their random nature is 
initially stemming from the microscopic random poten- 
tial V. 



Vst [F] oc exp 



and the Gaussian distribution 

^ I dy[dyFiyf^ (FIO) 

is a steady-state solution of (F8), i.e. dtPst [F] = 0, 
provided that A^^ = ^ (see Ref. 29). This condi- 
tion comes solely from the counterbalance of the dif- 
fusive term j^dyF{y) and the stochastic term y sF{y] ■ 
As for the contribution of the non-linear KPZ term 
— [dyF{y)]^, it disappears^'' under the boundary con- 
dition F'{y)\y±ac = 0. 

There is no 'time' dependence in the FP equation (Fl), 
but the Gaussian PDF (FIO) is a steady-state solution 
only at t = cx) for our physical definition of the pseudo 
free-energy (12) which satisfies 

dyFy{t,y) = dyFti,{t,y)+dyFy{t,y) = JL+r^y{t^y) 

(Fll) 

and it then becomes equivalent to the two normalized 
Gaussian PDF: 

m = ^exp(-^ / dy [dyF{y)f \ (F12) 



Ci(A) 
1 

C2(A) 



exp 



dy Hy)] 



(F13) 



with C2(A) defined by J 'Dri{y) ■ Vq [rj] ~ 1 and similarly 
for Ci(A). Choosing A^^ — ^ and the boundary condi- 
tion F'{y)\y±rx, = '']iy)\y±oo = 0- the FP equation even- 
tually yields for the Gaussian distributions: 

dtP'c [F] =-^0 [F] ■ log {C^ (A)7^g [F] } (F14) 



dyS{0)+\og{C,{X)P°,[v]} 



the equivalent equation for dtpQ [i]] being apparently ill- 
defined in the functional derivative framework, due to an 
additive divergent constant <5(0). So the Gaussian distri- 
butions (F12) and (F13) (with = 0) become a steady- 
state solution only at infinite 'time' in order to cancel the 
whole contribution (F14), since the value of A has been 
fixed as a 'time '-independent constant. 

Actually the chosen boundary condition is the only one 
physically possible at all 'times': 

dyFy{t,y)\y±oo = Vv{t,y)\y±oc = (F15) 

since by construction Vq [F] and Vq [?/] penalize the 
functions F{y) and r]{y) whose fiuctuations extend too 
much in the y-direction. The boundary condition at 
y — ±cx) is thus free from that point of view but should 
also be compatible with r]y(t, y) — (32). 



2. Steady-state solution at ^ = 



Steady-state solution of the linearized FP 
equation at ^ > 



For an uncorrelated disorder, the disorder correlator 
reduces to a normalized Dirac (5-function R(^=o{y) = 6{y) 



The Gaussian PDF (F12)-(F13) can be generalized 
with the introduction of the correlator R~^{t,y) whose 



38 



functional inverse is defined by (20): 

Vg [v,t] cx e-3/rfad!/''?fe)fl-Ht,la-!/'l)')(a') (F17) 

with the proper 'time'-dependent normalization of PDF 
such that / VF{y) ■ Pa [F,t] = J Vi^{y) ■ Vq [lA = 1. 

These expressions actually describe correctly the distri- 
butions corresponding to the linearized problem, namely 
to the equations (F1-F6) where the KPZ quadratic con- 
tributions to the functionals Q are set to 0. Indeed, the 
solution in the fields Fy, Fy and r]v of those equations is 
linear in the disorder potential V{t, y), whose distribution 
is Gaussian, implying that the distributions of the fields 
are themselves Gaussian. We refer the reader to Ap- 
pendix E for motivations to study the linearized prob- 
lem and for a solution leading to the two-point correlator 
-R''"(i, y) given in (40-41), and actually providing the full 
'time'-dependent distributions through (F16-F17) with 
R = R^"\ 

At ^ — the solution of the full problems (F1-F6) is 
not Gaussian at finite 'time', as known from the exact so- 
lutions^^'^'*''"'^'^^ but becomes Gaussian at infinite 'time', 
as detailed in the previous subsection. At ^ > this last 
property does however not hold anymore. Indeed, antic- 
ipating slightly on Sec. VB, if the steady distribution of 
r] were Gaussian then the three-point function ^3(^,0) 
would be zero by parity in the field 77, and from (91) the 
infinite- 'time' limit of the height R{t, 0) of the correlator 
would be the same in the linearized and in the original 
problem, which is not true. 

Another way to illustrate this fact is to try to solve the 
FP equation (F8) inserting the Gaussian Ansatz (F16). 
Using the fast decay of R{y) and R'{y) at large y for the 
vanishing of boundary terms, we obtain from the right 
hand side of (F8) 



dtPa [F,t] 
Vg \F,t] 



dydyF"iy)R-\t,\y-y\) 



xji^y dy'J^{t,\y - y'\)F"{y') - 2g [F;t,y]j (F18) 
with the definition: 



where M(t; j/', y) = M{t] y' — y) is a translation-invariant 
symmetric functional operator which combines R and R. 
Since (F21) is valid for any function F' decaying fast 
enough at infinity, solving this equation amounts to can- 
celing the operator M{t\ y' , y). After some manipulations 
aiming at casting the functional equation M(t;y',y) — 
into a diagonal form, the linearized form of the 'fiow' 
equation (34) on R{t, y) is precisely recovered, namely 

dtRit, y) =^dlR{t, y)~\ {Rit, y) + dy [yR{t, y)] } 



Dmy) 



(F22) 



In the process, the divergent part on the last line of (F18) 
was discarded. This fiow equation was obtained in Ap- 
pendix D using Ito's lemma, including the non-linearized 
case, and without having singular terms to discard (which 
we attribute to an artifact of functional calculus in the 
computation above). The infinite- 'time' steady-state so- 
lution thus verifies: ^dyR{oo,y) — DR'^{y) — which 

implies directly the expected result R{oo,y) = ^R^{y). 
It ensures that T{t,y) — >■ {^)~^6{y) in (F19) as t goes 
to infinity, which actually prevents the divergent term to 
appear in (F18). The finite-'time' solution is studied in 
Appendix E. 

Taking however the corresponding steady-state distri- 
bution (F16) with R{t,y) = ^R^{y) as trial steady so- 
lution for the full equation (F8) yields a remaining term, 
arising from the non-linearity, cubic in F (i.e. not of the 
form (F21)), which vanishes only at ^ — 0. 



Appendix G: Scaling laws in the convention of 
mathematicians 



The DES model defined in Sec. II A depends on the 
four independent parameters {c, D, T, with the elastic 
constant c being fixed independently from the temper- 
ature in the usual convention considered by physicists. 
However, an alternative convention used by mathemati- 
cians is c = T, whose consequences on the scaling argu- 
ments of Sec. IIIC and the saddle-point arguments of 
Sec. Ill D are discussed in this appendix. 

One is interested in the generic scaling of the prefactor 
of the roughness B{t;c,D,T,£^) in the random-manifold 



T{t,yi-y2)= / dysR^yi-ysDR {t,\y2-y3\) (F19) regime (of roughness exponent Ci 



RM 



-)• 



On the other hand, the left hand side of (F8) yields, 
differentiating with respect to 'time': 

ldydy'F'{y)d,R-\t,\y-y'\)F'iy') 

(F20) 

Considering first the linearized case, the identification 
of (F18) and (F20) yields, upon appropriate integrations 
by part, an equation of the form 



B{t; c, D, T, Anuic, D, T, 0^''^'' 



(Gl) 



The scaling in temperature of the prefactor is described 
by the thorn exponent ]d defined by A{c, D,T,£^) ~ T^^. 
We have derived in section IIIC from a scaling analy- 
sis that, depending on the temperature regime with re- 
spect to Tc = (^cD)^^^, the expressions of the prefactor 
Arm{c,D,T,£,) are 



dydy'F'{y')M{t;y',y)F'iy)^0 



(F21) 



(T»Te: 



D 



A 



RM 



(T<Te) 



(G2) 
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While lyj^^" = — I at high temperatures, the existence 
of the microscopic length ^ > alters the value of the 
thorn exponent to Ja^^"'"" = at low temperatures, even 
though this exponent describes large-scale properties of 
the polymer. 

In this appendix, we determine how those exponents 
change when taking the particular convention c = T, 
often chosen in the mathematics community - prompting 
us to denote by a subscript 'm' the observables defined 
with this convention e.g. 



(G3) 



Physically, the choice c — T amounts to render the elastic 
weight temperature-independent, (3) becoming 



w^^(ti,yi) 



a(ti)=ai 

Vye 

y(0)=0 



{dtv) 



+ T^V(t,y{t))^ 



(G4) 



Here ^ only tunes the relative importance of disorder 
with respect to elasticity. The parametrization c = T 
also arises in the continuum limit of the discrete simple 
SOS directed polymer model'^^'^^ and is thus of interest 
to analyze numerical results of this system. 

Before handling the different limits with respect to T 
of (G4) in a functional integral saddle-point approach 
similar to that of section HID, we first recall some re- 
sults on the asymptotics of integrals with one variable. 
The aim is to determine a (logarithmic) equivalent at 
large p of integrals of the form I{p) — J dy /(?/)e^''^^^^. 
The following result holds: if g{y) has a unique, finite, 
minimum value reached in y* , then 



Hp) " fiy*)e-P''^y*^ 



(G5) 



Here y* is the point (or one point) where the minimum of 
g{y) is reached, and is thus by definition independent of p. 
Powerlaw corrections in p may arise from the integration 
of fluctuations around y*, but they disappear e.g. in 
ratios of the following form (see (72) and (78) for the 
DP): 



I dy f{y)e-Psiy) fj^y-y-paiy') 
J dy e-P3(y) ^ 



-pa(y*) 



fiy*) (G6) 



The existence of the finite minimum is crucial, as illus- 
trated from the derivation of Stirling's formula for the 
equivalent of the factorial. Starting from 



g{y) = — logy reaches no finite minimum on R+. On 
this simple example, the clue is to rescale y hy a factor p 
{y — py) and write instead of (G7) 



p\ — p 



dy e'P^yP 



+1 



dy e 



-p{y-logy) 



now with f{y) = 1 and g{y) = y - 



- log y which is minimal 



in = 1 one obtains correctly''^ p\ ~ pP+^e"^. Note 
that, coming back to the initial variable y, we see that the 
optimum y* of g{y) in (G8) was not finite but diverging 
to infinity as y* = py* = p for p — >■ oo. In other words, 
the rescaling y = py in (G9) allows to find the optimal y 
at the correct scale in the large parameter p. 

Consider first the well-controlled high-temperature 
regime. Since the rescaling of section III C is at c = T = 
1, it is compatible with the mathematician's convention 
and one can export directly (54)-(55) imposing c = T 



' t ■ 1 1 ^ 



_D2 , lath 



D 



)(G10) 
(Gil) 



The regime f ^ 
with = 
limit and consists in replacing 
In the large-time limit {t ^ 



C^hiT), or equivalently T:^T^, 
describes the high-temperature 



by in (GIO). 
C(T)), this yields 
and the high- 



temperature thorn exponent is thus — " — |. 

The low-temperature regime is however less direct to 
handle, since the rescaling (56-57) of section III C is not 
at c — T and thus cannot be directly exported to the 
mathematician's convention. Anyway as first choice the 
rescaling a = ^, 6 = allows to rescale at ^ = 1 and to 
respect the mathematician's convention: the elastic term 
is unchanged (c ~ T) in the weight 



/■ti/5^ r( 



-\2 rpm 

+ ^Vi(t,j,(t)) 



W^ih,yi)'^ I Vye 

JyiO)=0 

(G12) 

where Vi{t,y(t)) = V{t,y{t))\^_-^ from which one 
reads 



dy e ^yP 



dy e ^e- 



y^piogy 



R+ 



(G7) 



one may be tempted to apply (G5) with 

f{y)^e-y g{y)^-\ogy (G8) 

which, assuming blindly that g{y) reaches a fi- 
nite minimum in y* would yield the wrong re- 
sult p\ e-y\P^°'iy\ The loophole here is that 



t . - 



1) 



(G13) 



However, the limit T — > cannot be taken by can- 
didly replacing by in (G13). This would lead to 
-^-2/3 g^jjj yield a corresponding zero thorn 



^RM I {t<s:t^ 

exponent, but this appears incorrect as we now discuss. 

T2 
T 

only in front of the disorder term, and not in front of both 



Indeed, the term ^ in the Hamiltonian in (G12) appears 
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contributions as in (72). In terms of path integrals 



ti ±1 

^) e Jo e 



JVy y{ 
. y(o)=o 



-[ '^dt^idtvf -'^l ^ dtVi{t,yit)) 
Jo g J Jo 



Vy e 

y(0)=0 



(G14) 



the large prefactor actually selects the path which 
minimizes the disorder contribution along the polymer 
trajectory, and not the full Hamiltonian as in (72), 
and this path is 'too anomalous' (no elastic constraint 
enforces it to stay in a bounded region as T — )■ 0). 
To contend with this singular limit, instead of start- 
ing with (G13), one may better work in the physicists 
convention starting from the scaling construction (56)- 
(57), where the low T behavior is controlled. One checks 
that the only possible rescaling of (G13) into a physicists 
roughness B{t; c, D, T, ^) satisfying c = D = ^ = 1 is: 



1/3 



C(r) 



D 



1/3 



(G15) 



B^{t-D,T,0 = B{t-T,D,T,0 



1,1, 



T 



1 (G16) 



This rescaling is similar in spirit to the rescaling y = py 
in (G9) for the saddle-point asymptotics study of pi : it 
allows to find the optimal path y{t) at the correct scale 

in the large parameter ^ encountered in (G13). The 
limit T in the mathematician's B"^{t; D,T,£^) coin- 
cides with the limit T — > in the physicist's roughness 

of (G16) since — - — r — -> 0. One reads from (G16) at 

(4T_D)3 

asymptotically large time, according to the known results 
Crm = 2/3: 



B^ihD,T,0 



C'^rm(i,1, 



(J TD)i/3 



,1 



t 1 ^^'^ 

(G17) 



In the low-temperature regime (T <^ T^) the physicist's 



A 



RM 



(l,l. 



r,l 



remains finite and goes to a T- 
0, as dis- 



T_D)i/3 

independent finite constant in the limit T 
cussed previously in Sec. HID. This finally yields 



(G18) 



At low temperature, the mathematician's thorn exponent 
]d is thus non-zero: t>RMl(T<T"')= ~ |- 



Appendix H: Set of parameters for the numerical 
simulations 



We have listed thereafter the parameters of the nu- 
merical study presented throughout Sec. IV, and initially 
defined in Sec. IV A. 

The discretization grids for the generation of the ran- 
dom potential and for the recording of the data are fixed 
once and for all to the following values for all the com- 
putations (the last three parameters are respectively the 
number of points for the linear grid in t, the logarithmic 
grid in t and the linear grid in y as defined in Fig. 4): 



£)giid 




fgrid 


nbptlint 


nbptlogt 


nbptliny 


4 


1 


2 


80 


80 


100 



The following sets of data have been generated, with 
the elastic constant being fixed at c = 1 and at different 
temperatures T. The number of configurations 'NconfV' 
per data set have fluctuated as a compromise between 
the convergence of the disorder average and a reasonable 
computation time, as already mentioned in Sec. IV A. 
The first column gathers the sets used for the study of 
R{t,y), Dt and £^t, whereas the second column corre- 
sponds to the additional sets which explore larger 'times' 
t < tm and have thus been used for the study of the 
roughness B{t): 



no 


T 




2/m 


NconfV 




Urn 


NconfV 


1 


0.35 


40 


50 


572 








2 


0.4 


40 


50 


1092 








3 


0.5 


40 


50 


728 








4 


0.55 


40 


50 


540 








5 


0.6 


40 


50 


654 








6 


0.7 


40 


50 


1281 








7 


0.8 


40 


50 


2665 








8 


0.9 


40 


50 


900 








9 


1 


40 


50 


1350 








10 


1.1 


40 


50 


1040 


100 


70 


518 


11 


1.35 


40 


50 


1050 


300 


60 


390 


12 


1.5 


40 


50 


1435 


200 


80 


1600 


13 


1.6 


40 


50 


1600 








14 


1.8 


40 


50 


1000 


800 


160 


159 


15 


2 


40 


50 


1750 








16 


2.25 


40 


50 


1300 








17 


2.5 


40 


50 


1900 








18 


2.75 


40 


50 


1000 








19 


3 


40 


50 


2300 








20 


3.2 


40 


50 


1000 








21 


3.5 


40 


50 


1825 








22 


4 


40 


50 


1100 








23 


6 


40 


50 


900 


80 


50 


280 
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Appendix I: Effective correlator of a 2D 
cubic-splined microscopic disorder 

As described in Sec. IV, a given configuration of the 
microscopic disorder V{t, y) is generated first by pick- 
ing up a set of random numbers {V, } on a grid of spac- 
ing l^f"^"^! Cf"^"^! with a normal distribution of variance 

D^"'^, and then by interpolating between the grid points 
with a 2D cubic spline (cf. Fig. 22 top). 

In the reduced case of a ID cubic spline between a 
set of points at a fixed position t on the grid, it is pos- 
sible to determine analytically the effective correlator of 
the ID cubic-splined random potential with a normalized 
function i?^"''''^^ (y) (cf. (1)). For the 2D sphne we then 
assume that the effective two-point disorder correlator is 
given by the translational-invariant: 



V{t,y)ViO,0)^D-R 



CubicS 



it) 



^CubicS (y) 



(11) 



with ^( = and £,y ~ by construction, and be- 
cause of the passage from the discretized yj to the con- 
tinuous variable y an amplitude D = £)gi'id^s"'^^gnd^ 

In practice, for 2n -I- 1 points indexed by j = —n, . . . , n 
on a grid of spacing ^, a random value Vj is attached 
on each site yj = of the grid following the statistical 



distribution V-j 



and VjV-j' 



-'jj' ■ 



A cubic spline of 



{Vj}_n<j<n ^ function V{y) which is a cubic polyno- 
mial on each lattice segment y S [yj, yj+i], continuous on 
each lattice site yj and with its first and second deriva- 
tives also continuous at y = yj . Combining the equations 
of the cubic-splined parameters for a given set {Vj} and 
the disorder average over those possible sets, we obtain 
the following symmetric effective correlator in the limit 
of n — >■ oo: 



R' 



CubicS 



(0 < y < 



7i iy- 



e 







{4-3V3)y^+^y + e 



<y < 



y s vj+i) 

J {y ~ yj) (y - yj+i) (Vs 
y-yi-i - y^{y~yj) 



(12) 



(13) 



where there is a distinction between the central segment 
5: 2/ < which contains the auto-correlation 



nCubicS 



(y = yj=()) 



^0 



1 



and the other segments yj < y < yj+i {j — ^ ■ ■ - n) with 
oscillations constrained by the cancellation at y = yj^j 



\y^y,^o)^VoV,^o^O 



Note that this correlator has been obtained exactly by 
averaging over disorder the cubic splines, whose coeffi- 
cients depend linearly in the random potential and thus 
allow an analytical computation of this average. 




y(o,j/)y(o,o) 




V{0,y)V{0,0) 



20 y 



SincG (+1.0) 
- CubicS (+0.5) 
8 Measured 

y 



FIG. 22. (Color online) Top: Visualization of a sin- 
gle disorder configuration V{t,y); the mesh accounts for the 
grid {yj}, and fluctuates smoothly. Bottom: Graph of the 
numerically-measured effective-ID cubic-splined disorder cor- 
relator (blue), with the 'SincG' and 'CubicS' (red, respectively 
dashed and continuous) superimposed with a slight transla- 
tion for more visibility. Inset: Full effective disorder correla- 
tor at flxed t — 0, computed over lO'OOO disorder configura- 



grid 



tions with D^"^ = 4, ^f"" = 2 and C. 



For the microscopic disorder correlator (II) the func- 
tion R^ is thus exactly known, but it is not the case 
for the effective disorder rj{t,y)rj{Q,Q) — D ■ R^{y). As 
explained in Sec. IV we try nevertheless to extract the 

effective parameters at fixed time |_D,f| with three dif- 
ferent fitting functions normalized to 1: a Gaussian (88) 
('Gauss'), a sinus cardinal with a Gaussian envelope (89) 
('SincG') and the exact cubic-spline correlator (I2)-(I3) 
('CubicS'). The comparison between those three options 
accounts for a consistency check for the determination of 
the effective parameters via the stability of the function. 
Having first tested this procedure on V{t,y)V{0,0) (II), 
by fitting it separately at fixed time t — versus at fixed 
position y = 0: 



Vit, 0)1/(0, 0) = [D- (0)] • (t) (X % it) 
V{0,y)V{0,0) = [D ■ R^, (0)] • R^^ (y) ex R^^ (y) 



(14) 



we summarize thereafter the values obtained for the am- 
plitude D and the correlation lengths ^t,y with respect 
to the grid parameters (those are kept at fixed value 
throughout all our computations). The errors are pro- 
vided by the nonlinear-regression procedure on a corre- 
lator obtained over lO'OOO disorder configurations: 
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D 




6 


i?i(0) 


Grid 


4 


2 


1 




Gauss 


10.08 ±0.08 


0.625 ±0.005 


0.312 ±0.003 


1 


SincG 


8.04 ±0.02 


2.020 ± 0.004 


1.012 ±0.002 


1 

Erf(7r) 


CubicS 


7.95 ±0.02 


2.006 ± 0.004 


1.005 ±0.002 


1 



As we can check both in this table and visually in Fig. 22 
(bottom), the sinus cardinal ijSmcG matches the main phe- 
nomenological features of the exact /j;CubicS^ ^ ^^le am- 
plitude of the central peak and the position of the first 
oscillation, and keeps satisfied the relations £,t y — Cfj)*^ 

and D = i:)g"d^grid^grid^ 



As for the Gaussian fit, it overestimates the amplitude 
of the peak and underestimates its typical variance. For 
the latter it is simply due to the geometrical definitions 
of the parameter ^ proper to each of these functions, 
and which can be merged up to a numerical constant 
(^CubicS/^Gauss ^ 3 ^ov the first, it stems from the 
fact that a normalized Gaussian function is used to fit 
a normalized function which has negative contributions: 
the integral of the central peak has thus an unnormal- 
ized integral larger than 1, and fitting it by a Gaus- 
sian naturally yields a larger corresponding amplitude 
pCubicS/£)Gauss _ q gg-j^ rpj^-^ discrepancy will actually 

be present in all the similar fitting procedure in Sec. IV. 
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